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Abstract 

We construct supergravity plus branes solutions, which we argue to be related to 4d A/" = 1 
SQCD with a quartic superpotential. The geometries depend on the ratio Nj/Nc which 
can be kept of order one, present a good singularity at the origin and are weakly curved 
elsewhere. We support our field theory interpretation by studying a variety of features like 
i?-symmetry breaking, instantons, Seiberg duality, Wilson loops and pair creation, running 
of couplings and domain walls. 

In a second part of this paper, we address a different problem: the analysis of the interesting 
physics of different members of a family of supergravity solutions dual to (unflavored) M = 1 
SYM plus some UV completion. 
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1 Introduction 

The AdS/CFT conjecture originally proposed by Maldacena pi and refined in |2l EI is one 
of the most powerful analytic tools for studying strong coupling effects in gauge theories. 
There are many examples that go beyond the initially conjectured duality and first steps 
in generalizing it to non-conformal models were taken in |3]. Later, very interesting de- 
velopments led to the construction of the gauge-string duality in phenomenologically more 
relevant theories i.e. minimally or non-supersymmetric gauge theories [Sj. 

Conceptually, the clearer setup for duals to less symmetric theories is obtained by breaking 
conformality and (partially) supersymmetry by deforming J\f = 4 SYM with relevant oper- 
ators or VEV's. The models put forward in [Sj and (even though there are some important 
technical differences) the Klebanov-Tseytlin and Klebanov-Strassler model(s) [Hj belong to 
this class. 

On the other hand, a different set of models, that are less conventional regarding the UV 
completion of the field theory have been developed. The idea here is to start from a set of 
Dp-branes (usually with p > 3), that wrap a g-dimensional compact manifold in a way such 
that two conditions are satisfied: first, one imposes that the low energy description of the 
system is (p — g) -dimensional, that is, the size of the g- manifold is small and is not observable 
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at low energies. Secondly, one also requires that a minimal amount of SUSY is preserved, in 
order to have technical control over the theory (for example, the resolution of the Einstein 
eqs. is eased). According to intuition, in this second class of phenomenologically interesting 
dualities, the UV completion of the (usually four-dimensional) field theory of interest is a 
higher dimensional field (or string) theory. There are several models that belong to this latter 
class. In this paper we will concentrate on the model dual to TV = 1 SYM jTj , that builds on 
a geometry originally found in 4-d gauged supergravity in [S]. It must be noted that all of 
the models in this category (and also those in the class described in the previous paragraph 
jSlini) are afflicted by the fact that they are not dual to the "pure" field theory of interest, 
but instead, the field theory degrees of freedom are entangled with the KK modes (on the 
g-manifold or with the modes in the latest steps in the cascade) in a way that depends on 
the energy scale of the field theory. The KK modes enter the theory at an energy scale which 
is inversely proportional to the size of the g-manifold and the main problem is that this size 
is comparable to the scale at which one wants to study non-perturbative phenomena such 
as confinement, spontaneous breaking of chiral symmetry, etc. Nevertheless, this limitation 
can be seen as an artifact of the supergravity approximation and will hopefully be avoided 
once the formulation of the string sigma model on these RR backgrounds becomes available. 
Many articles have studied different aspects of these models. Instead of revisiting the main 
results here, we refer the interested reader to the review articles 0. 

An important characteristic of the different models mentioned above, is that these super- 
gravity backgrounds are conjectured to be dual to field theories with adjoint matter (also 
some quiver gauge theories are in this class). A problem that remains basically unsolved 
is how to deal with field theories containing matter fields transforming in the fundamental 
representation. It is a very natural problem to tackle if one is interested in making contact 
with phenomenological theories, like QCD. 

There were some attempts to study this problem. To begin with, the "quenched" approx- 
imation (when the number of flavors is negligible compared to the number of colors, hence 
using a probe brane approximation in the dual picture) was studied in detail and in the 
different models, leading to a set of interesting results, regarding breaking of symmetries, 
meson spectrum and interactions. This line was developed in detail; for a list of references, 
see the citations to the paper that initiated this approach [TUj . 

The need to go beyond this "quenched" approximation motivated some authors to search 
for more complicated backgrounds including the effects of a large number of flavors jTTj. 
These new backgrounds, that typically depend on two coordinates, are afflicted (in the 
case of four-dimensional fleld theories), by the usual problems caused by the addition of D7- 
branes: a codimension-two object with an associated conical singularity, a maximum number 
of D7-branes (or O7 planes) that can be added in a compact space, etc. There is another line 
of research that used non-critical string theory ^21 1^ to construct duals to interesting fleld 
theories. The use of non-critical strings cures the problem of the existence of extra massive 
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states (the KK modes pointed out above), but it typically relies on backgrounds which have 
large curvature, with the consequent problem that solutions need to be stringy-corrected 
(see P5 for some work in this direction). The results obtained are nevertheless indicating 
that this approach has some potential. 

Perhaps being conservative, one could say that the results regarding the duality with 
flavored fleld theories using the backgrounds of [HllISl, are not as spectacular and clear as 
the ones obtained with duals to fleld theories with adjoint flelds only. All this indicates that 
a different approach is necessary, which might try to combine the successes of the previous 
ones and to avoid their inconveniences as much as possible. 

In this paper we propose a way of flnding a dual to A/" = 1 SQCD-like theories using 
critical strings, focusing our attention here on type IIB string theory. The basic idea will be 
to add flavors to the background in [ZI, by using different technical points developed in [12] 
and in non-critical string approaches ^21 d] ■ Below, we describe our approach and give a 
guide to read this paper. 

A second part of the paper, namely section |H1 deals with backgrounds dual to different 
UV completions of minimally SUSY Super- Yang-Mills. For these geometries, unlike in the 
model [Zj, the dilaton does not diverge in the UV. We study relevant aspects of the dual 
gauge theory. 

1.1 General Idea 

Let us describe here the procedure we will adopt to add a large number of flavors Nf to 
a given supergravity background, more particularly to the one constructed with wrapped 
D5-branes and conjectured to be dual to A/" = 1 SYM plus massive adjoint matter [Zj. Even 
though we will concentrate on this particular case, the procedure described below could be 
used to add many flavors in different backgrounds dual to A/" = 1 , A/" = 2 SYM in 2 + 1 and 
3 + 1 dimensions, etc. 

We want to preserve SUSY in order to solve BPS eqs rather than Einstein eqs; since the 
original background preserves four supercharges, we cannot break more SUSY in duals to 
four-dimensional fleld theories. It is then necessary to flnd a non-compact and holomorphic 
two-cycle (S2), where we can place Nf D5-branes that share the 3+1 gauge theory dimensions 
with those Nc D5-branes that generated the original background [Zj. These surfaces must be 
holomorphic to preserve SUSY, and non-compact so that the symmetry added by the flavor 
branes is global (in other words, the coupling of the effective four-dimensional fleld theory 
on the flavor branes gl = gl/Vol(T,2), vanishes) [T7j . 

The problem of flnding non-compact SUSY preserving two-cycles in the geometry of [7] 
was solved in We will use a particular solution found there in order to place many 

D5-branes. This new stack of Nf branes is heavy, so it will backreact and the original 
background will be modifled. We can think of an action describing the dynamics of the 
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backreacted system that reads ^ 

•S" = Sub + S flavor (1-1) 

Basically, we have added the open string sector in S flavor- The procedure might be thought 
of as consisting of two steps; first we take the background [7], where the Nc color branes 
have been replaced by a flux, that we represent by Sub in On top of this, we add a 

bunch of Nf flavor branes and find a new background solving new BPS (and Einstein) eqs. 
that encode the presence of the flavor branes. The part of the action corresponding to 
the flavors will be the Dirac+WZ action of the Nf five-branes. In Einstein frame it reads 

Sfiavor = T,y2(- I d'xe^.r^) + [ P[Ce]] , (1.2) 

where the integrals are taken over the six-dimensional worldvolume of the flavor branes Ai^ 
and ^(6) stands for the determinant of the pull-back of the metric on such worldvolume. We 
will take these D5's extended along six coordinates that we call Xq, Xi, X2, x^, r, -0 at constant 
6,(p,6,(p. These supersymmetric embeddings were called cylinder solutions in |15j . 

The action Sub is ten-dimensional in contrast with S flavor, that is six-dimensional since 
the "flavor branes" are localized in the directions 6, 6, ip, ip. So, trying to find solutions to 
the equations of motion derived from ()l.ip . will involve writing Einstein eqs with Dirac delta 
functions on the coordinates where the branes are localized. Thus, the solution will depend 
on {6, 9, ip, if) and in this respect, this solution would be similar to those found in [IT]. This is 
a very hard problem in principle and it is nice to notice that the (six- dimensional) non-critical 
string approach of ^2] to A/" = 1 flavored solutions is not afflicted by this technical difficulty, 
since S flavor is basically a cosmological term ^. We will circumvent this difficult technical 
obstacle by following a procedure first proposed in ^Hl in the context of flavor branes. In 
this paper, an eight-dimensional non-critical string action was used and an ingenious trick 
developed. Indeed, it was found in jTHI that one could solve Einstein eqs without Dirac 
delta functions if one smears the very many Nf branes in the extra directions (in our case 
0,0, if, (p). 

Since Nf and N^ are very large numbers of the same order, interchanging the sum in 
fll.2|l by an integral over the extra directions will produce a fully ten-dimensional action, 
erasing the explicit dependence on the extra coordinates. This gives rise to some new global 
symmetries. We will comment on the interpretation of this procedure from the dual field 
theory perspective in section 15.11 

-'^It is known that it is not possible to write a polynomial action for IIB supergravity due to the self-duality 
condition, we nevertheless will deal with solutions that have F^, = 0, so, what we have in mind in this case 

is an action of the form Sub ~ J d}'^X\/9{R — \id4>Y — ■ 

^But as we mentioned, the non-critical string approach is seriously afflicted by string corrections to the 
gravity approximation; nevertheless, the AdS solutions in (12j are likely to persist after these corrections. 
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We proceed in the rest of the paper by first describing in detail the background in jTj and 
a singular version of it (that is easier to deal with and illustrates the technical details). We 
will then propose a new deformed background, where the fingerprint of the fiavor branes will 
be the explicit violation of the Bianchi identity (a sort of Dirac-string like singularity). We 
will find BPS eqs describing this situation and in appendix |X1 will develop a superpotential 
approach and compare with a purely IIB SUSY approach to this problem. The outcome 
being that the BPS eqs obtained from the superpotential coincide with those coming from 
imposing vanishing of the variations of the gravitino and dilatino, Stp^, 6X, in ten dimensions. 

Then, we apply a similar approach to the non-singular case, obtain a set of BPS eqs, find 
some interesting asymptotic solutions with numerical interpolation and study their (strongly- 
coupled) dual gauge theory predictions. 

The content of the paper is the following: in section 2, we present the dual to A/" = 1 SYM 
plus adjoint matter that will be the arena on which we will construct fiavored solutions. The 
presentation is detailed enough to make the paper self-contained and the reader familiar 
with these results may skip it. In section IHl we will provide a derivation of the fiavored- 
BPS equations for the singular background, that solve the Einstein eqs derived from 
This derivation is complemented in appendix using a superpotential approach. The 
presentation in this section will be very detailed, and it was written in order to describe 
and explain in a simpler context the procedure we will follow in the physically interesting 
non-abelian case. 

In section m we will write BPS equations for the fiavoring of the non-singular background 
and, as an interesting particular case, we will also deal with the case Nf = 2Nc that presents 
unique features. Then, we study in detail the asymptotics of the solutions and provide a 
careful numerical treatment to the eqs derived in this section. In section El we first provide 
arguments to describe the field theory dual to our backgrounds, that indicate that (at low 
energies) we are dealing with Af = 1 SQCD plus a quartic superpotential for the quark 
superfield. Then, we will initiate the study of the field theory properties of these fiavored 
solutions, most notably, we will analyze Wilson loops and SQCD-string breaking, instanton 
action, theta angle, beta function, f/(l)_R symmetry breaking, domain walls and Seiberg 
duality. 

In section IHl we will come back to the interesting case Nf = 2Nc. Here, we will analyze 
distinctive gauge theory features; most notably (non)-confinement (screening of quarks), 
U{1)r symmetry preservation and Seiberg duality, that becomes particularly interesting in 
this case. We will also comment on finite temperature aspects of this field theory. 

In section a different approach to flavored backgrounds will be introduced. We argue 
that it might be possible to account for the violation of the Bianchi identity induced by the 
addition of the flavor branes, by turning on non-trivial fluxes which are not present in the 
unflavored solution of 

In the second part of the paper (as a byproduct of our results above), in section |H1 we 
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will consider the particular case Nf = (which in the following we will call "unflavored 
case" ) and find a one parameter family of deformed (non-singular) solutions that correspond 
to different UV completions of A/" = 1 SYM. We present a gauge theory interpretation of 
this family of solutions and study many field theory aspects as seen from it; most notably, 
confinement, k-string tensions, PP-waves, rotating strings, and beta function, pointing out 
in each case the similarities and differences between different members of the family and the 
solution in [7]. 

Section El is left for conclusions and possible future directions. With the aim of making 
the main text more readable, we wrote many appendixes, where we have relegated lots of 
technical details. 

Reader's guide Given that this is a considerably long paper, but that different parts may 
be read almost independently, we feel it is useful to write a guide in order to help the reader 
find the results of her/his interest. Readers interested in the construction of the fiavored 
backreacted solutions can concentrate on sections El and El supplemented with appendices 1X1 
and IHl and also read section |3 Readers interested in the field theory features reproduced 
by the fiavored solutions can go directly to sections El and El and just look at the explicit 
expressions for the backgrounds (sections 14. 2114. 5|) when needed. Finally, who is interested 
in the deformed unfiavored solutions can read, in a self-contained way, sections 14.11 and El 
supplemented with appendices |Bl and |El 

2 The Dual to A/" = 1 SYM 

We work with the model presented in (the solution was first found in a 4d context in 
[S]). Let us briefiy describe the main points of this supergravity dual to A/" = 1 SYM and 
its UV completion. We start with A^^^ D5-branes. The field theory that lives on them is 6D 
SYM with 16 supercharges. Then, suppose that we wrap two directions of the D5-branes 
on a curved two-manifold that can be chosen to be a sphere. In order to preserve some 
fraction of SUSY, a twisting procedure has to be implemented JH] and actually, there are 
two ways of doing it. The one we will be interested in in this paper deals with a twisting 
that preserves four supercharges. In this case the two-cycle mentioned above lives inside a 
Calabi-Yau 3-fold. The corresponding supergravity solution can be argued to be dual to a 
four- dimensional field theory only for low energies (small values of the radial coordinate). 
Indeed, at high energies, the modes of the gauge theory explore also the two-cycle and as the 
energy is increased further, the theory first becomes six- dimensional and then, blowing-up 
of the dilaton forces one to S-dualize. Therefore the UV completion of the model is given by 
a little string theory. 

The supergravity solution that interests us, preserves four supercharges and has the topol- 
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ogy M}'^ X M X S"^ X S"^. There is a fibration of the two spheres in such a way that A/" = 1 
super symmetry is preserved. By going near r = it can be seen that the topology is M^'^ x 5^. 
The full solution and Killing spinors are written in detail in ^Sl- Let us write some aspects 
for future reference and to set conventions. The metric in the Einstein frame reads, 

dslo = a'gsN.e^ \ —^dxl ^ + dr^ + e^^ ( dO^ + sin' Od^^ ) + - {ui - A'f 1 , (2.1) 
L a'QsNc ' ' 4 J 

where is the dilaton. The angles 9 G [0, vr] and ip G [0, 27r) parametrize a two-sphere. This 
sphere is fibered in the ten-dimensional metric by the one-forms {i = 1,2,3). They are 
given in terms of a function a{r) and the angles {6, (p) as follows: 

= -a{r)de, A^ = a{r)smedip, A^ = -cosOdip. (2.2) 

The a)j one-forms are defined as 

cDi = cos ipdO + sin sin Odp , 
U2 = — sin ipdO + cos ip sin 6d(p , 

(Da = dip + cosOdif . (2.3) 

The geometry in ()2.ip preserves four supercharges and is non-singular when the functions 
a(r), h{r) and the dilaton (f){r) are: 

2r 

a[r) - 



e'^ = r coth 2r 



sinh 2r 

r2 1 



sinh' 2r 4 ' 



^-24> ^ ^ (2.4) 

sinh 2r 

where 0o is the value of the dilaton at r = 0. The solution of type IIB supergravity includes 
a RR three-form F(3) that is given by 

-^F(3) = -\{u,-A')A{u,-A')A{us-A') + ^ F« A - A'^ ) , (2.5) 

where F° is the field strength of the SU{2) gauge field A"", defined as: 

= dA^ + ^eabcA'AA' . (2.6) 
The different components of read 

F^ = -a'drAd9, = a sinOdr A dp , = (1 - ) smOdO A dip , (2.7) 
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where the prime denotes derivative with respect to r. Since (i-F(3) = 0, one can represent F(3) 
in terms of a two- form potential C(2) as F(3) = dC(2)- Actually, it is not difficult to verify 
that C{2) can be taken as: 



(2) 



ip ( sin 6d6 A d(f — sin 6d6 Adip) — cos 6 cos 6d(f A dip 
a{d9 AuJi — sin 9dip A 0J2 



{21 



The equation of motion of -F(3) in the Einstein frame is d{^e'^*F(^)^ = 0, where * denotes 
Hodge duality. Let us stress that the configuration presented above is non-singular. 

For future reference, let us quote here the asymptotic expansions of the functions a(r), 

g2h(r) g20{r) ^^^-^ ^ ^ g 



a(r) ~ 1 r -|- 



1 + 



and for large values of the radial coordinate. 



,2r 



air) ~ 4re"^'' + ... 



_ g20o. 



4v^ 



(2.9) 



(2.10) 



The BPS eqs that the configuration in ()2.1|l - ()2.8j) solves (for a complete derivation from the 
spinor transformations of IIB sugra, see appendix A in J3]) read 



h' 



1 

Q 
1 

2Q 
2a 



-2h 



„2h 



16 



(a^ - 1)^ 



+ 1 + 



ia' - 1? 



e^h + ^ (a^ - 1) 
4 ^ ^ 



with 



The Einstein eqs read 



Q 



+ i ^2h (^2 + 1) + ^ (a2 _ ly 
2 16 



(2.11) 



(2.12) 



The Maxwell equation as quoted above reads 6*^*^(3) 
dF3 = 0. Both of them are solved by (EH) -(ESI). 



and the Bianchi identity is 
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Finally, in the next section, as a warm up example, we will add many flavor branes to a 
particular singular solution of the system ()2.1H) . which is characterized by a(r) = 0, e^^ = 



r, e 



2<f>-2<f>o 



and a RR potential 



ip ( sin 9d9 A dip — sin 6d6 A dip) — cos 6 cos 6dip A dip 



(2.14) 



(2.15) 



Notice that at r = this background presents a (bad) singularity that is solved by the 
turning-on of the function a{r) in the full solution (|2.1|) - (j2.8|) . On the field theory side, this 
way of resolving the singularity boils into the phenomena of confinement and breaking of 
the R-symmetry. Since they will be useful in the remainder of this paper, let us summarize 
some aspects of the dual field theory. 



2.1 Dual Field Theory 

In jTj the solution presented in ()2.1|) - ()2.8|) was argued to be dual to A/" = 1 SYM plus some 
KK massive adjoint matter. The 4D field theory is obtained by reduction of D5-branes 
on S"^ with a twist that we explain below. Therefore, as the energy scale of the 4D field 
theory becomes comparable to the inverse volume of S"^, the KK modes begin to enter the 
spectrum. 

To analyze the spectrum in more detail, we briefly review the twisting procedure. In 
order to have a supersymmetric theory on a curved manifold like the S'^ here, one needs 
globally defined spinors. In our case the argument goes as follows. As D5-branes wrap 
the two-sphere, the Lorentz symmetry along the branes decomposes as 5*0(1,3) x S0{2). 
There is also an SU{2)l x SU{2)ji symmetry that rotates the transverse coordinates and 
corresponds to the R-symmetry of the supercharges of the field theory on the D5-branes. 
One can properly define A/" = 1 supersymmetry on the curved space that is obtained by 
wrapping the D5-branes on the two-cycle, by identifying a U{1) subgroup of either SU{2)l 
or SU{2)r R-symmetry with the 5*0(2) of the two-sphere. To fix notations, let us choose 
the U{1) in SU{2)l. Having done the identification with 5*0(2) of the sphere, we denote 
this twisted f/(l) as U{1)t. 

After this twisting procedure is performed, the fields in the theory are labeled by the 
quantum numbers of 5*0(1, 3) x U{1)t x SU{2)fi. The bosonic fields are (the a indicates an 
adjoint index) 

= (4,0,1), <|.'^ = (1,±,1), r=(l,±,2). (2.16) 
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Respectively they are the gluon, two massive scalars that are coming from the reduction of 
the original 6D gauge field on S*^ (explicitly from the and Ag components) and finally 
four other massive scalars (that originally represented the positions of the D5-branes in the 
transverse W^). As a general rule, all the fields that transform under U{1)t, the second entry 
in the above charge designation, are massive. For the fermions one has, 

A'^ = (2, 0, 1), (2, 0, 1), = (2, ++, 1), (2, — , 1), r = (2, +, 2), (2, -, 2). (2.17) 

These fields are the gluino plus some massive fermions whose U{1)t quantum number is 
non-zero. The KK modes in the 4D theory are obtained by the harmonic decomposition 
of the massive modes, ^, and ip that are shown above. Their mass is of the order of 
^KK = (^0/52)"^ oc -^-^jjf- A very important point to notice here is that these KK modes 
are charged under U{1)t x U{1)r where the second U{1) is a subgroup of the SU{2)r that 
is left untouched in the twisting procedure. On the other hand, the gluonic gauge field and 
the gluino are not charged under either of the f/(l)'s. 

The dynamics of these KK modes mixes with the dynamics of confinement in this model 
because the strong coupling scale of the theory is of the order of the KK mass. One way 
to evade the mixing problem would be to work instead with the full string solution, namely 
the world-sheet sigma model on this background (or in the S-dual NS5 background) which 
would give us control over the duality to all orders in a', hence we would be able to decouple 
the dynamics of KK-modes from the gauge dynamics. This direction is unfortunately not 
(yet) available. Meanwhile, in [TH] a procedure was developed to determine when a field 
theory observable computed from the supergravity solution is affected by the presence of the 
massive KK modes or is purely an effect of the massless fields A^,\. 

In order to get a better intuition of the dynamics, one can schematically write a lagrangian 
for these fields as follows: 

L = -Tr[^Fl^ + i\D\ - (D^$,)' - {D^i.f + m{iD - M^k)^ + Mlj,{il + $2) + V[i, $, ^ 

(2.18) 

The potential typically contains the scalar potential for the bosons, Yukawa type interactions 
and more. This expression is schematic because of (at least) two reasons. First of all, the 
potential presumably contains very complicated interactions involving the KK and massless 
fields. Secondly, there is mixing between the infinite tower of spherical harmonics that are 
obtained by reduction on 5*^ and S*^, see pUl 1^ for a careful treatment. 

3 Adding Flavor Branes to the Singular Solution 

In this section we will add flavors to the particular solution ()2.14|) - ()2.15p . Even though there 
is little physical significance for this solution, the point of this section is to illustrate in 
detail the type of formalism we use and its remarkable consistency. Readers more interested 
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in physically more relevant aspects of this paper should perhaps skip this section in a first 
reading, but since the formalism and technical subtleties are quite involved, we decided to 
spell them out explicitly here in a simpler context. 



3.1 Deforming the Space 

The first step will be to deform the background ()2.14j) - ()2.15j) . First, we will study the 
deformation of the backgrounds without the addition of flavors and then we will treat a 
deformation due to the presence of the Nj flavor branes. For this we propose a set of 
vielbeins given by, 

e^* = eW-, e'^ = e^dr, e*^ = e^+'^dO, e'^ = e^+'^sin 
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e = , e = , e = . 3.1 

2 2 2 ^ ^ 

that leads to the metric (Einstein frame): 

2g{r) 2k{r) 

ds" = e^f^'\dxl^+dr^+e^''^'-\de^+siYi^ ed^^) + '^^{ul+ul) + '-^{u:^^^^ . (3.2) 

Notice that compared to ()2.14|) . we took a'gs = 1, while Nc has been absorbed in e^^, e^^, 
e^'^. In our configuration, we will also have a dilaton and a RR three- form 

0(r), F(3) = -2iV,e-3/-29-fcei A A + ^e'^Z-^'^-^e^ A A e^ (3.3) 

where N^, the number of color D5-branes, comes from the quantization condition: 

/ F(3) = N^T, . (3.4) 

The 5*^ on which we integrate is parameterized by 6,ip,il), so only the first term in ()3.3p 
contributes. 

We plug this configuration into the IIB SUSY transformations (see eq flB.3|) ). The projec- 
tions for the Killing spinor are: 

Trme = e, r^,0<^3e = e, e = ie* . (3.5) 

(The sign of the two first expressions can be changed, but we fix it in order to match the 
full non-abelian case). After some algebra, the following eqs. arise, 

4/ = (3.6) 

h' = + ^e-2^+' (3.7) 

g' = + (3.8) 

k' = h^^e-^^-^ -N,e~^'J-^ -^e-^^+^ -e-^3+^ + 2e-^ (3.9) 

(j)' = —Nce-^''-'' + Nce~^3-'' (3.10) 
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3.2 Flavoring the Space 



We want to backreact the space with D5-flavor branes. Let us try to do this following 
the procedure introduced in jT2j, i.e adding an open string sector to the gravity action. It 
is important to remark that the construction of involved non-critical strings and thus 
curvatures of order of the string scale. This fact hindered the possibility of finding reliable 
quantitative results from a supergravity approximation. On the other hand, we are using 
ten-dimensional string theory and weak curvature can be obtained by taking gsN^. ^ 1 as 
usual. The action of the system is: 

•S* = Sgrav + S flavor (^-H) 

where Sgrav (Einstein frame) is given by: 



9 - " / d^^T.r^ 

'-'grav o 2 / u/ j-y y 



i(9,0)(9'^0)-le<^Ff3) 



(3.12) 



whereas S flavor is the Dirac+WZ action for the Nf D5-flavor branes (we take the worldvolume 
gauge field to be zero). In Einstein frame: 

Sfiavor=T,J](- [ (fxe^.r^)+ I P[C^]] (3.13) 

where the integrals are taken over the six-dimensional worldvolume of the flavor branes Aie 
and ^(6) stands for the determinant of the pull-back of the metric in such worldvolume. We 
will take these D5's extended along xq, ■ ■ ■ X3,r,ip at constant 6,ip,6,(f. Notice that this 
configuration preserves the ?7(l)i? symmetry associated to shifts in ip. This is one of the 
embeddings which were called cylinder solutions in jTH]. Moreover, notice that this brane 
configuration makes clear the need of the deformed ansatz ()3.2j) where k{r) ^ g{r) (compare 

to dHD). 

Following a procedure analogous to |16j we think of the Nf —>■ oo branes as being homo- 
geneously smeared along the two transverse parameterized by 6, (f and 6, (p. We will 
elaborate on how the smearing may affect the dual field theory in section 15.11 The smearing 
erases the dependence on the angular coordinates and makes it possible to consider an ansatz 
with functions only depending on r, enormously simplifying computations. We use the same 
ansatz ()3.2|1 for the metric and also fix = 4/. We have: 

-T^J^J^ d'xe'^^f^) ^ j Asin^sine~e^y^ (3.14) 

Nf 

^[^e] ^ jt^J Vol{y,)AC(^e) (3.15) 



13 



where we have defined ^^0/(3^4) = sin 9 sin 9d9 A dip A d9 A d(f and the new integrals span the 
full space-time. We will need the expressions (we choose a' = gs = I)- 

= fr, ' TT^ = /o \7 (3.16) 

(27r)5 2k^^o) (27r)7 

Let us turn our attention to the effect of the WZ term of the flavor brane action ()3.15|) . 
Since it does not depend on the metric nor on the dilaton, it does not enter the Einstein 
equations. However, it alters the equation of motion for the 6- form C(6), which before was 
d * F(7) = (i-F(3) = 0. Now we have a source term so^: 

ciF(3) = ^NfVol{y4) = ^Nf sin 9 sin 9d9 Ad^Ad9 Ad^ (3.17) 

Notice that this particular form for the violation of the Bianchi identity is an effect of the 
smearing: we replace the sum of delta functions on the position of each of the Nj branes by 
a constant density, so there is a continuous distribution of charge which acts as a source for 
the RR form. We can slightly modify ()3.3p to solve ()3.17p "^: 

JSl ~ ~ Nf - N 
F(3) = — ^ sin 9d9 Ad(p A {dip + cos 9d(f) ^— — - sin 9d9 A dip A {dip + cos 9d(p) (3.18) 

On the other hand, we still have (iF(7) = 0. In order to flnd the new system of BPS 
equations, we recompute ()B.3j) using the modified RR three-form fleld strength ()3.18|1 . It is 
quite straightforward in this way to get: 

h' = ^{N,- Nf)e-^^-'' + (3.19) 

k' = h^N,-Nf)e~^^''^ -N,e~^3-^ -^e-^^'+^ -e-^3+k j^2e-^ (3.21) 
0' = ~{N^- N})e~^^-^ + N^e-^3-k (3_22) 

Nevertheless, it is important to point out that, apart from the IIB sugra action, we also have 
the flavor branes action, so, in order to preserve supersymmetry, it is necessary that their 



•^In general, if for a form — c?A(„_i) there is an action — 2^ / ^/WW^ + / G A A, the equation of 
motion for the form reads d * F — sign{g){—l)^^"^^G. In this case, the relevant part of the action (go to 
string frame for this computation) reads: — ^J^^ J ■y/jgjFjy^ + J ^0/(3^4) A C(6) so the equation of 



motion is ^^l^ d * f (y) — ~j^^Vol{y4). Taking into account i^(3) = — *F(7) and (|3.1t)|l we arrive at H3.17|l . 
''More generally, we could consider: 

Nc + N' ~ ~ Nf - Nc - N' 
F(3) = — - sin 0d0 A dip A {dip + cos Odip) sin OdO A dip A {dijj + cos 9d(p) 

but N[. is just a shift in Nc- As explained below (|3.4|l . it is the coefficient of the first term the one identified 
with the number of colors. 
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action preserves /t-symmetry. As stated above, this is indeed the case in this construction. 
In appendix we reobtain this system using the so-called "superpotential" approach and in 
section ing we check that this first order system automatically implies that the set of second 
order equations of motion are satisfied. 

The system ()3.19p - ()3.22j) has a very simple solution when Nf = 2Nc that we present in 
section 14.51 and discuss in section El In this work, we will not undertake the study of the 
system for Nf ^ 2Nc. 

Using this toy-example, we have clarified our proposal. In order to add backreacting 
flavors, we will consider a deformed background, solution of the eqs of motion derived from 
the sugra action plus an action like ()3.14|) - ()3.15|) . These last terms come from the action 
of the (susy preserving) flavor branes. We modify the RR form so that the failure of the 
Bianchi identity indicates the presence of the smeared flavor branes. After that, we impose 
vanishing of the IIB SUSY variations. This will produce a set of BPS eqs that satisfy the 
Einstein and Maxwell eqs and ensure the susy of the whole construction. 



3.3 The Einstein Equations 

In order to check for good the consistency of the approach, in this section we will check that 
the solutions found above indeed satisfy the full set of equations of motion. 

First of all, we have already guaranteed the equation for the RR-form (see ()3.17p and 
(I3.18|l ). The equations for the worldvolume fields on the flavor branes are also satisfied since 
K-symmetry is preserved. We still have to check the equations for the dilaton and the metric. 
The relevant action is given by the sum of ()3.12|) and ()3.14|) . since ()3.15|) does not depend 
on the dilaton nor on the metric. The equation for the dilaton is: 

-7^9, {g^''^^)dA) - ^e'^Ff - ^e^^^^ sin0sine~ = (3.23) 
V-5'(io) J-2 8 V-fl'(io) 

whereas the equations coming from variations of the metric are (compare with (j2.13p ): 

(3.24) 

where the energy-momentum tensor generated by the flavor branes comes from the variation 
of the lagrangian of the flavor branes ()3.14j) : 



avor 



^flavor = - sin 6* sin 6*6 2 -y/ -^(6) (3.25) 



and therefore is: 



_ 2/t(xo) SC flavor _ Nf n-n'^i f M g^/^ V ~^(6) 

flavor - , ^ - — smysmy-e2 o^Of^g,^^. 
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where span the brane world volume. In components, after lowering the indices: 



rpflavor 






rpflavor 

rr 


^f ~2h-2g 
2 


rpflavor 
^ipip ~ 


^f 2k-2h-2g 
8 




rpflavor 
tpip ~ 


_^^2k-2h-2g Q 
8 


rpflavor 

W ~ 


_^^2k-2h-2g ^^^2 
8 


0, 


rpflavor 


_^^2k-2h-2g Q 
8 


rpflavor 
(p(p — 


8 


0, 


rpflavor 


_^^2k-2h-2g^^^Q^^^~Q 
8 



Now it is straightforward to check (using Mathematica) that the equations (jH.19j) - (j3.22j) 
imply that ()3.23p and ()3.24p are satisfied. 

To summarize, we have explicitly showed in this section that a consistent procedure to 
add flavor branes to a given background is to consider the BPS eqs obtained by imposing the 
vanishing of SUSY transformations for the gravitino and dilatino in a system of a "deformed 
spacetime" (hke ()3.2j) with respect to ()2.14|l ) and modified RR forms so that the failure of 
the Bianchi identity is indicating the presence of the smeared flavor branes. In the next 
section, we will apply this proposal to add flavor branes to the background dual to A/" = 1 

SYM (ini)-(EIHl)- 



4 Adding Flavor Branes to the Non-singular Solution 

In this section we will solve the important problem of adding Nf flavor branes to the back- 
ground dual to N = \ SYM ()2.H) - ()2.8|) ^ . As pointed out in the introduction, the resolution 
of this problem is instrumental in the duality between string theory and M = 1 SQCD-like 
theories. We will be more sketchy here than in sectionlHland just present the ansatze and cor- 
responding BPS equations for different cases. Many technical details are left for appendix iBl 
Our procedure (as explained in detail in section El), will be first to propose a deformation 
of the background dual to A/" = 1 SYM ((2111) -(IZHl)- Then, we will obtain the BPS eqs that 
describe the deformation due to the presence of flavor branes, that will be extended along 
the directions Xq, Xi, X2, X3, r, ip and smeared over the directions {6, 9, ip, (p). As before, these 
flavor branes will be sources for RR forms that will flux the deformed background and the 
mark of their existence as singular sources will be the violation of the Bianchi identity. 

4.1 Deforming the Space 

We start by proposing a deformation to the background ()2.1|l - fl2.8j) . The ansatz we use below 
is a subcase (with no fractional branes if(3) = F(5) = 0) of that first proposed in j2SI and 

^In the probe approximation also the paper |22j discussed the problem. 
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further analyzed in j2l|. We borrow some results and notation from j24j. Although it is 
not the main motivation of this work, this unfiavored setup encodes some interesting physics 
which we will analyze in section |H| We consider the Einstein frame metric: 



+ 



e2/(r) 

4 



dxl^ + dr' + e^'^'^de' + sin' Od^') + 
[{uji + a{r)def + {Cj2 - a{r) smedipf) + 



,2fc(r) 



{uj3 + cos9dipf . (4.1) 



The vielbein we consider for this metric is the straightforward generalization of ()3.1|) by the 
inclusion of the a{r) dependence in and e^. Apart from the dilaton 

= 4/, (4.2) 

we also excite the RR 3- form field strength which we take to be of the same form as ()2.5|) : 



^(3) - -J- 



(oji + h{r)d9) A {uo2 - h{r) sinOdif) A (cDs + cosOdip) + 



b'dr A {-de A cJi + sin Odip A UJ2) + (1 - h{rf) sin OdO A dip A 



(4.3) 



The condition (iF(3) = is automatically ensured by this ansatz. It is useful to rewrite this 
expression in terms of the vielbein forms: 



N, 



\3) 



2N,e-''^-^'-'e' Ae^Ae' + tl^h' e'^f-^-^e' A {-e' A + A e^) + 



^_3/-2h-fc(^2 _ + ly /\ /\ gS ^ N,e-''f-^'-3-^{h - a) (-e^ A + A e'^) A e^(4.4) 



We now analyze the dilatino and gravitino transformations ()B.3jl (in this case, the super- 
potential approach is far more complicated because, as we will see, there are algebraic con- 
straints). After very lengthy algebra (that is explicitly reported in appendix Ej) 7 we obtain 
a system of BPS eqs and constraints, some of which can be solved, leaving us with two 
differential eqs for a and k. We define: 



-^^'^dr = dp . 



(4.5) 



The differential equations read: 



dpk 



2p coth2p 



e^'' (acosh2p - 1) 



Nr 



sinh 2p 



+ a{2p — a sinh 2p) 



2(l + a2-2acosh2p)-i 
{-1 + 2pcoth2p) 



— a sinh 2p(a cosh 2p — 1) + (2p — 4ap cosh 2p + -— sinh 4p) 

(4.6) 
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The equation for / = | reads 

Qr_ (-1 + acosh2p)^sinh"^(2p)(-4p + sinh4p) 

~ 4(l + a2-2acosh2p)(-l + 2pcoth2p) ' ^ ' 

and can be integrated once a(p) is known. The functions b{p), g{p), h{p) can be solved to 
be: 



b(n) = 2,^^ ^COsh(2p)-l 

sinh(2p)' ^acosh(2p)-l 



g2h_ ^(2acosh(2p)-l-a2). {U 



2g 



4.2 Flavoring the Space 

We consider the same ansatz for the metric: 



dx'l^ + dr' + e^^^'\dd'' + sin" Od^') + 



edipf) + {Cj^ + cos edipf\ . (4.9) 

and the dilaton is still (j4.2p . 

One should now think of introducing backreacting flavor branes along the lines of sectional 
We incorporate smeared flavored D5-branes on the non-singular solution and the analysis 
of the equations is quite similar. We again consider D5's extended along xq, . . . , 0:3, r, ■?/', 
each at constant 6', ip, 6, (p. The analysis in ^H] shows that these branes preserve the same 
supersymmetry as the background for any value of the angles 6', yj, 6', One can therefore 
think of smearing along this space^. As in section the Bianchi identity gets modifled to 
jHZ!): 

rfF(3) = ^NjVoliy^ = ^Nf sin 6 sin Ode A d<f A dO A d<f (4.10) 
We solve this by adding to the 3- form written in ()4.3|) the same as in ()3.18|) : 

pflavor _ _ ^ Q^Q ^^^^ ^ PQg ^^n^ 

Written in flat indices, the 3-form now reads: 

F(3) = -2N^ e-^^-^3-k^i A A + ^ t'e'^^-'^-^e'' A (-e^ Ae^ + e^ A e^) + 
^ ^ ^-3f-2h-k ^^2 _2ab + l-x)e' Ae^ Ae' + 

+ N, e-'^f-^-^-\h - a) (-e^ Ae^ + e^ A e^) A (4.12) 



^Notice that in the unflavored construction the spheres parameterized by 9, if, 0, (p do not play an impor- 
tant role from the point of view of the IR A/" = 1 SYM theory. Since we want to add flavor to this theory, we 
can expect that the effect of the smearing will not alter dramatically the resulting IR A/" = 1 SQCD theory. 
We will comment more about this in section II 
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where we have defined 

x^f (4,13) 



We now insert the expression ()4.12|1 for F(3) in the transformation of the spinors equations 
(jB.3jl . We find a set of BPS eqs: seven first order equations and two algebraic constraints 
(for details, the reader is referred to appendix |Bj. This system'' can be partially solved, 
leaving us with (here we use the definition ()4.5j) ): 

^ _ (2 - x) p _ 2b cosh 2p - 2 + X 

sinh(2p) ' 2 acosh2p— 1 

^2h ^ ^(2acosh(2p) - 1 - a^) (4.14) 

plus two coupled first order eqs for a(p), k{p) and a differential equation for / = | that can 
be integrated in terms of a(p). These equations read {k = k — \ log Nc): 

2 / 2e2^ -a;(acosh2p- 1)2 o ■ . ...rN 

d.a = ^ ^ + smh 2p - 2ap (4.15) 

^ 2pcoth2p-ll 2-x sinh2p h hj k j 

2 f 2e^^ + X 

dok = -. r- asinh2p(acosh2p — 1) + 

(2pcoth2p-l)(l-2acosh2p + a2) I 2-x ^ ' 

2sinh 2p 



and 



+2p(a^ cosh 2p - 2a cosh 2p + 1) (4.16) 

2p 



4p + sinh4p + 



(-1 + acosh2p) sinh"2(2p) 

4(1 + a2 - 2a cosh 2p) (-1 + 2p coth 2p) 

4 -1 

+4ap cosh 2p - 2a sinh 2p - -a(sinh 2p)^ (4.17) 

(2-x) J 

The equations ()4.14|1 - ()4.17|1 constitute an important result of this paper. In order to sum- 
marize the flavored setup, the metric is given by ()4.9|) . the dilaton by ()4.2j) and the RR 
three-form is the sum of ()4.3|) and ()4.11|) (or, in flat indices, it is ()4.12|) ). while all the func- 
tions in the ansatz are determined by ()4.14j) - ()4.17j) . One can check (using Mathematica) that 
these conditions solve the full system of second order equations given by ()3.23|) . ()3.24|) and 
d * F(3) = 0. Notice that the expressions ()3.27p are still valid in this case. So, if we can flnd 
solutions to ()4.15|1 - ()4.17|1 . we will have the string dual to a family of SQCD-like theories with 
Nc colors and Nf flavors; we now turn into this. 

^These expressions are clearly not valid when Nf — 2Nc. We will deal with this special case in section l^T^ 
and appendix IdI 
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4.3 Asymptotic Solutions {Nj ^ 2Nc) 



Unfortunately, it was not possible for us to find exact explicit solutions to the system ()4.15p - 
()4.17|1 . We study the problem in the following way, we first find asymptotic solutions near 
p — s> oo (the UV of the dual theory) and near p = (the IR of the dual gauge theory). Then, 
we study a numerical interpolation showing that both expansions can be joined smoothly; 
for many purposes, this is as good as finding an exact solution. 



4.3.1 Expansions for p ^ oo 

Let us disregard exponentially suppressed terms, and assume that a = e~^''j where j is some 
function that can be Laurent expanded in p~^. Then, when x < 2, the large p expansion 
reads: 

a = e-'" ((4 - 2x)p + | + ||^p-^ + 0{p-^)^ +... (4.18) 

When a; = we have a = Ape"'^^ + . . . which agrees with the usual solution (|2.1|) - (|2.8|) . see 
()2.10|) for comparisons. The other functions read: 



dj = --—p-^ + ... (4.19) 
Equation ()4.19|) is clearly not valid for x > 2. In this case we find a different expansion: 

The other possible behavior at p — cxd is a geometry asymptoting to Minkowski times the 
deformed conifold. Since we will not use such solutions in our physical interpretation, we 
relegate their description to appendix O 
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4.3.2 Expanding around p — > 

A series expansion can be found for tlie functions a{p),k{p) that solve the BPS eqs near 
p = 0. This is actually a two-parameter family of solutions labelled by two free numbers 
which we denote ci and C2- The solution for the functions reads in this case 

. ^ 2 s (80 - 40a; + 9a;c?) . 



12(2 -a;) 



2k_.r I 2 ^ 3ciC2a; 3 C2(256 - 256x + x^(64 + 27c^)) , 

e — iVc I C2P + —, -7-p H — -— r;^ p + 



2(x-2)'^ 48(2 -x)2 

2„ /2(2-x)l a; 8(2 - a;) 

\ 6Ci p 2 9ci 

2k „/'2(2-iO, I, 8(2-1), 



,2</.-2</.o 



Some comments about this expansion are in order. First, we should notice that ()4.2ip above 
does not reduce to (j2.9|) when x = 0. This is not very surprising, since the behavior of 
gauge theories in the IR changes radically in the presence of massless flavors. Second and 
perhaps more important, we notice that this solution is singular. Indeed, when computing 
the Ricci scalar when p ^ 0, one gets R ~ p~^ and there does not seem to exist a choice of 
the constants (ci,C2) that avoids this. 

The presence of a singularity might be source of concern and cast doubts on the reliability 
of the solutions. An important point that we would like to emphasize, though, is that there 
is an important and significant difference between the singularity in 1)4.211) and others present 
in the literature, as the one in ()2.14j) for example. There is one criterium developed in 
(see also [211 other criteria) that suggests a way of deciding when an (IR) singularity in 
a supergravity solution should be accepted or rejected as unphysical. This criterium is quite 
easy to apply and coincides with other different criteria developed in previous literature. It 
consists in analyzing the (Einstein-frame) component of the metric Qu- If this is bounded, 
the singularity should be accepted as a good background. The idea in this criteria is that an 
excitation travelling towards the origin (that is a low energy object on the gauge theory dual) 
will have an energy as measured by an inertial gauge theory observer given by E = \/\gtt\Eo 
(with Eq the proper energy of the object). So, if gu diverges (the singularity is repulsive), 
the object gains unbounded energy from the field theory inertial observer point of view and 
this is to be considered as unphysical. 

According to this criterium, the singularity encoded in the solution 1)4.21)) above should 
be accepted as being physically relevant; indeed gu = e"^/^ = e'^°^'^{l + ...), see ()4.21)) . This 
means that even though the (good) singularity we find might be a signal of some IR dynamics 
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that we are overlooking, we can still perform computations to try to answer non-perturbative 
field theory questions. As we will see in section El the background matches several gauge 
theory expectations even when we are performing computations near p = 0. 

To complete the study of these solutions, it is worthwhile to do a detailed numerical 
analysis to which we turn now. 



4.4 Numerical Study {Nf < 2Nc) 

In order to get a qualitative idea of the behavior of the solutions, we have studied numerically 
the equations ()4.14|) - ()4.17p with the initial conditions ()4.21|) near p = 0. We also impose the 
large p asymptotic behavior ()4.19|) . what gives a relation between the parameters Ci and C2- 
We also fix C2 = 9cf/16. This amounts to requiring that the domain wall tension scales with 
N(. — ^ (see section 13^ . However, the qualitative behavior of the plots does not depend on 
this constraint. With these conditions, all the functions are uniquely determined for a given 
value of X = up to the constant 0o which we will set to zero. Some results are reported 
in figure m 

Notice that the a; ^ limit is not continuous (at least in the small p region) as it is 
apparent from the expansion ()4.21|) . In order to show it graphically, we have included in 
figure H a plot of the the function a near p = 0. 



4.5 The Simple Nj = 2Nc Solution 



From the expressions ()4.14|1 - ()4.17|1 . one can readily see that Nf = 2Nc is a special point and, 
in fact, the manipulations leading to some of those expressions are ill-defined in this case. 
Nevertheless, from the first equation in ()4.14|1 we expect to have 6 = 0. Then, considering 
the system (jRTT^ - fjRTTjl . it is easy to see that it can be solved taking a = b = and ^ = 1, 
B = (however, this is not the only possibility, see appendix |D)) . In fact, a = b = leads 
us back to the construction of section IH.2I For the particular Nf = 2Nc case, there is an 
extremely simple solution, in which the other functions satisfy: 

e'' = ^, ^'' = r=^' ^'' = ^- = 4/ = 0o + r (4.22) 

where we have introduced the constant < ^ < 4. The resulting Einstein frame metric, 
dilaton and RR three-form read: 



3 + Nc(dr^ + ^{de'^ + sin 



+ 



4-e 



{de^ + sin^ ed<p 



+ - {dip + COS 9d(f + cos OdifY 



\3) 



2 



sin 9d9 Adip + sin 9d9 Ad<^) A{dip + cos 9dip + cos 9d(p) (4.23) 
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Figure 1: Some functions of the flavored solutions for Nf = 0.5Nc (thin solid lines), Nf = 
l.2Nc (dashed lines), Nf = l.GNc (dotted lines). For comparison, we also plot (thick solid 
lines) the usual unfiavored solution. The graphs are e"^^, e^^, e"^^, 0, a and a zoom of the plot 
for a near p = 0. In the figures for a we have also plotted with a light solid line a = — tttt- 



We will discuss the gauge theory associated to this background in section El 
4.5.1 A Flavored Black Hole 

In section 14.3.21 we discussed the fact that the < Nf < 2Nc solutions present a curvature 
singularity, which is good according to the criterium of having bounded gu- The Nf = 2Nc 
case presents the same kind of singularity at r — oo. However, the simplicity of this 
metric allows us to check that this singularity also satisfies an alternative criterium, that 
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basically states that a gravity singularity is good if it can be smoothed out by turning on 
some temperature in the dual field theory and thus finding a black hole solution [221 ■ So, in 
order to construct the black hole, let us first make a change of variable, 

r = 2logz, (4.24) 

and define: 

jr=i_(^E^y (4.25) 
Keeping the same dilaton and three-from but deforming the metric to: 

/ 4 1 

- J^dt^ + dxl + dxj + dxl + NA—r~^dz^ + -{dO'^ + sin^ Qd^"^) + 

-XdP + sin^edif"^) + i(# + cosedif + cos edifY)] (4.26) 



7 2 ^ 

ds = e ^ z 



provides a regular solution of the Einstein equations ()3.23|) . (j3.24p . As a technical comment, 
notice that, due to the change in the metric, the expressions for T^""^"^ and Tj^°-^°'^ get 
multiplied, with respect to those in (j3.27|) by JF and respectively. The rest of (j3.27j) 
remains unchanged. 

Let us stress that this solution has its own interest since it is, to our knowledge, the first 
solution dual to a field theory in four dimensions with adjoints and fundamentals at finite 
temperature*. We will explore some of the gauge theory aspects of this solution ()4.26p in 
section [HI 

Finally, to close this section, we mention a numerical solution described in detail in ap- 
pendix O This solution for the case Nf = 2Nc, appears after a careful analysis of the BPS 
eqs and has remarkable features. We comment more on it at the end of section IHl 



5 Gauge Theory Aspects and Predictions of the Solu- 
tion {Nf ^ 2Nc) 

In this section we will study the family of gauge theories dual to the solution(s) we have 
presented in the previous section. We will first try to give arguments towards a definite 
lagrangian dual to the background (|4.9|) - (|4.17|) and then we will study different strongly- 
coupled aspects of these field theories to support our proposal. These aspects include R- 
symmetry breaking, pair creation and Wilson loops, instantons, Seiberg duality, domain 
walls and asymptotic beta functions. 

^In an interesting work |27j, the addition of finite temperature to a three dimensional gauge theory with 
adjoint and fundamental fields was studied. 
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5.1 General Aspects of the Gauge Theory 



Below we will present some arguments to motivate a lagrangian description for the field 
theory dual to our backgrounds ()4.9j) - ()4.17|l . We will argue that the dual field theory is 
related in the IR to = 1 SQCD plus a quartic superpotential for the quark superfields. 
Some aspects related to the smearing of the flavor branes and to the six- dimensional little 
string UV completion of the theory are not completely under control, and in particular we 
could expect that some operator we are overlooking might still be present and slightly deform 
the IR dynamics of the theory. Nonetheless we are confident that the result to which we 
arrive is robust and we present many tests in the following, which support our interpretation. 

Our argument starts by reminding the discussion in section 12.11 that the field theory dual 
to the unflavored solution of jZj is six-dimensional M = 1 SU (Nc) SYM compactified on a 
twisted two-sphere in such a way that the low-energy theory is four-dimensional and A/" = 1 
supersymmetric. The spectrum of the corresponding U{1) theory is related to that of the 
four- dimensional U (M) Af = 1* theory, that is A/" = 4 SYM with equal masses for the adjoint 
chiral multiplets, in its Higgs vacuum 20j. Let us briefly review the main features of this 
four- dimensional theory before arguing how to extend this remarkable result to the large 
case. The superpotential of the TV = 1* theory can be written (in A/" = 1 notations) as 

and therefore the F-term equations of motion read 

[$i,<l>j] = i/ieyfc$fc (5.2) 

It is immediate to recognize an 5^(2) algebra in this equation. This allows for non-trivial 
solutions to ()5.2|) where the adjoint fields $j are taken to be proportional to the generators 
jf"'* of any irreducible representation of su{2). For any dimension n there is only one such 
representation and therefore the classical vacua of this theory correspond to the partition of 
M into positive integer numbers [2H] 

M 

^nkn = M (5.3) 

n=l 

where fc„ is the number of times the n-dimensional irreducible representation appears in the 
expectation value of the adjoint fields. Any specific such vacuum has residual gauge group 
Yin U{kn)- To obtain a vacuum with gauge group f/(l), one has to take /cm = 1 and all other 
knS to zero, whereas to obtain a vacuum which has exactly SU{Nc) residual gauge group, 
one needs to start from SU{M) A/" = 1* with M an integer multiple of Nc, M = mNc, and 
to take km = Nc in ()5.3p . The expectation values of the adjoint fields in this Higgs vacuum 
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satisfy the relation 

^1 + ^1 + ^1 = ^2 (5.4) 

that is they describe a (fuzzy) two-sphere. 

In the particular case U{1), it was shown in (201 ^^^^ matching of the spectrum to 
that of the compactified six-dimensional theory we started with, is exact in the M — oo 
limit. Because of the ways the U{1) and the SU{Nc) vacua are built, it is natural to argue 
that the results of [201 HI] extend to the more general SU{Nc) case, that is the spectrum 
and lagrangian of the km = Nc Higgs vacuum of A/" = 1* is exactly the same as that of the 
six- dimensional SU{Nc) SYM wrapped on a twisted two-sphere when m —>■ oo. 

Let us consider now the case we are interested in, that is the addition of flavors to the 
compactified six-dimensional theory. When we add Nf flavor branes into the background, 
we are effectively adding a set of massless chiral multiplets transforming in the fundamental 
and anti-fundamental of the SU{Nc) gauge group. Let us call these multiplets Q = {S, q) 
and Q = {S, q). These fields are massless because the flavor branes are at zero distance from 
the color branes that generated the background. 

A generic lagrangian for this new system could then be written as 

j d^e (^$t^^e^$^^ + Q^e^Q + Q^e'^Q) + j d?e [W^ + W) ] (5.5) 

where W is a chiral superpotential describing the mass terms and self-interaction of the 
adjoint Kaluza-Klein states, and their coupling to the fundamental fields. Because of the 
argument above, we can think of the Kaluza Klein states as the three massive adjoint fields 
of the M = 1* theory. The most natural way to couple the fundamental fields to an adjoint 
state is mediated from the only allowed term in = 2 SQCD, which is the one typically 
appearing in intersecting branes setups. It reads 

W = kQ^Q (5.6) 

with K a coupling constant whose value will not affect our discussion. In a theory like the 
A/" = 1* in a Higgs vacuum described above, if the fundamentals couple to a single adjoint 
field through ()5.6|) . the SU{2) global symmetry that rotates the three adjoints breaks down 
to a U{1) subgroup. The term ()5.6|1 breaks explicitly also the SU{Nf) x SU{Nf) flavor group 
to its diagonal SU{Nf) subgroup. 

The explicit flavor group breaking is exactly matched in our dual string picture by the 
presence of flavor D5-branes alone, rather than of flavor D5 and anti-D5 branes. This 
feature is in contrast with other models of theories with fundamental degrees of freedom 
where the chiral symmetry breaking is spontaneous and flavors are added to the picture via 
both D-branes and anti D-branes, which reconnect in the IR to account for the spontaneous 
breaking [29j. Since in our case the breaking of the global symmetry is explicit rather than 
spontaneous, we do not expect to find Goldstone bosons in the spectrum. 
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Regarding the breaking of the SU{2) to f/(l), a httle more care is required since the 
smearing of the flavor branes restores the SU{2) symmetry in our model. Let us start, 
therefore, from a different distribution of the flavor branes to try to understand what is 
happening. We could smear the flavor branes on the {6, ip) two-sphere, while putting all of 
them on a single point on the {9, (p) directions, for simplicity let us say at 6* = 0. Since the 
exact unflavored solution of j7j is invariant under an SU{2) acting on the ip, 9 and p angles, 
this new configuration breaks this SU{2) background isometry to U{1). We claim that this is 
the configuration dual to the fundamental- adjoint coupling From the field theory point 

of view this would probably be the most natural brane configuration to consider, but as we 
already mentioned above, finding a solution to the Einstein and flux equations for this less 
symmetric configuration would probably turn out to be technically impossible. This is the 
reason why we considered a smeared configuration for the flavor branes. Let us try, therefore, 
to understand what this means in field theory terms. We know from ()5.4|1 that the three 
massive adjoint fields describe a two-sphere (in the limit M ^ oo it is a smooth sphere). 
Therefore, we can take them to be parameterized by two angles {9 and p), by the usual 
projections over the coordinate planes of of a two-sphere of radius = ^^{M"^ — l)/4 
centered at the origin. We will denote these adjoints parametrizing the S'^{9,(p) as "^l^"*^)- 

Then, we propose that the fundamental-adjoint coupling corresponding to our model is, 
schematically: 

W ~ . d9dp.r.~9Q + $1^-),^ + $1^-),^) Q (5.7) 

Let us quickly go through now the low-energy effective lagrangian corresponding to the 
superpotential ()5.7|) . In general, if the parameter /i is much bigger than the scale of strong 
coupling Ksqcd we can integrate out the KK fields from ()5.5|1 to end up with an A/" = 1 
SQCD-like theory that looks 
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j d''9 (gte^Q + Qte^Q j + J £9 + W] (5. 



where the effective superpotential W" will depend on the way the adjoints couple to the 
fundamentals and to themselves in ()5.5|) . For our case ()5.7|) . after integrating out the massive 
KK fields, we might expect to obtain some combination of the quartic operators ()5.6p : 

W ^ -{TriQQQQ) + (TrQQf) (5.9) 

and possibly some more complicated operator due to the cubic coupling of the adjoint KK 
fields. 

Could we take the KK masses to be very large (// oo), we would end up with pure 
SQCD. Unfortunately, there is no known way to sensibly separate fi from Aggcd without 
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requiring the full string dynamics^. In fact, one can expect that taking /i — oo requires 
reducing the "external" space (that is topologically 5*^ x S^) to a stringy scale. In this 
context, it is very suggestive that our construction has some similarity with the non-critical 
string model of Klebanov and Maldacena p2j, which, in fact, should be the limit of our 
construction when continuously increasing /i — >■ oo. Let us emphasize once again that, in the 
process, the curvature of space-time reaches the string scale and therefore one cannot expect 
to get reliable results from supergravity (although it has been shown that it is possible to 
get some qualitative information from that kind of setups using a gravity action as a toy 
model [lailSllini)- 

We can now complete the comparison of the isometrics of the field theory and of our 
background. We have already analyzed the matching of the flavor symmetry and the SU{2) 
Kaluza-Klein symmetry. The field theory at weak coupling also has a U{1)b x ^(l)i? symme- 
try, where U{1)b is an exact baryonic symmetry, while U{1)r is the anomalous i?-symmetry. 
The unbroken baryonic symmetry could be associated with rotations along whereas the 
anomalous i?-symmetry corresponds to shifts of the coordinate ip- As we will analyze in 
detail in section EiH the supergravity background reproduces faithfully the breaking pattern 

Finally, let us mention that in the special case A^/ = 2Nc we can reasonably claim that 
all the extra uncontrolled terms appearing along with ()5.9p in the effective superpotential 
W' can be turned off, giving on the string side our simple solution ()4.23|1 . As we show in 
section this background has very special features which we relate to the scale invariance 
of the dual field theory. On the other hand, for N j < 2Nc, we could not find a way on 
the gravity side to turn off these extra terms, as the comparison between the solutions of 
section 14.41 and of appendix O to the simple solution ()4.23|) seems to confirm. Nonetheless, 
as all the rest of this section shows, the string dual perspective of many strong-coupling 
effects, remarkably agrees with expectations from the field theory described by ()5.8j) - ()5.9|l . 
and allows therefore to make interesting predictions for its non-perturbative dynamics. 

As a concluding remark, it would certainly prove very interesting to try to match with our 
approach, the results of the careful analysis done in on the vacuum structure of A/" = 2 
SQCD with supersymmetry broken to A/" = 1 by a mass term for the adjoint superfield, 
which is closely related to the theory ()5.8|1 . 

^This is the typical problem that afflicts all the supergravity duals to confining field theories. This will be 
cleanly solved when a worldsheet CFT is found for these models lifting the limitation a' 0. Meanwhile, 
from a supergravity perspective, the KK modes can be disentangled form the gauge theory dynamics following 

[HEniini- 
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5.2 Wilson Loop and Pair Creation 

One natural question is what happens to the gauge theory Wilson loop, now that we have 
massless flavors. In principle, as in QCD we should not observe an area law, but the SQCD- 
string should elongate until its tension is equal to the mass of the lightest meson, and then 
break. So, if we find that the (very massive) quarks that we introduce in the system feel each 
other up to a maximal distance only, this will be indication that we are seeing a phenomenon 
like pair creation^". In order to study this, we will follow a very careful treatment to compute 
the Wilson loop in gravity duals developed in (for a good summary of the results, see 
pages 19-25 in ^^). As usual, we propose that the Wilson loop for two non-dynamical quarks 
(strings stretching up to p — oo) separated a distance L in the gauge theory coordinates, 
should be computed as the action of a fundamental string that is parametrized by t = r, 
Xi = 0", p = p(cr) [34]. By solving the Nambu-Goto action, one obtains a one-parameter 
family of solutions depending on an integration constant which we will define to be po (the 
minimal p reached by the string). We convert to string frame the metric fl4.9p and use ()4.2p 
and ()4.5j) . Then, the length and energy (renormalized by subtracting the infinite masses of 
the non-dynamical quarks) read: 



where pi is a cutoff related to the mass of the heavy quarks that can be taken to infinity 
smoothly. These integrals can be performed numerically. In order to do so, we have to fix 
the parameters ci and C2 of the expansion ()4.2H1 in order to give initial conditions for the 
equations (j4.15|) . (|4.16j) . Imposing the large p asymptotic behavior (|4.19|) gives a relation 
between them. Moreover, we also fix C2 = 9c^/16 so that the domain wall tension scales 
with Nc — ^ (see section . However, this condition is not essential for the qualitative 
behavior. The numerical result is reported in Figure El There is a finite value of po, say 
Pq, for which the solution reaches a maximum Lmax (at this point, ^ = ^ = 0). There 
is also a minimum value of the length = t^VOsNcQ' . Our interpretation is that Lmm is 
related to the KK scale at which the UV completion sets in, and, therefore, the background 
loses predictivity. For a range of Lmin < L < Lmax, there are two solutions of the Nambu- 
Goto action. The physical one (the one of lower energy), corresponds to po > Pg which in 
particular implies that the string never approaches p = and therefore it does not approach 
the curvature singularity. 

As explained above, when one stretches a flux tube in a confining theory with dynamical 
quarks, the tube must break at some finite length when there is enough energy to create 

^°In [33], screening due to backreacting flavors was found as a change in the sfope of Regge trajectories. 




(5.10) 
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Figure 2: Energy vs. length of the Wilson loop for jf- = 0.2 and = 0.5. Energies 
are in units of '^'^"^f ^'^ , whereas lengths are in units of yfgjT^ .There is a value of the 
integration constant po for which the string reaches its maximum L and E. Longer strings 
are not solutions of the Nambu-Goto action. We interpret it as string breaking due to 
quark- ant iquark pair creation. 



a light meson which screens the color charge. The geometry can reproduce this behavior 
(the breaking happens at L^ax) because the backreaction of the flavors has been taken 
into account. The graphs qualitatively match those found on the lattice in similar (non- 
supersymmetric) setups jSZ]. Thus, our solution is showing that there is pair creation and 
giving a prediction for the lightest meson in terms of the parameters of the model. The 
fact that Ljjiin and L^ax are roughly of the same order reflects the known fact that, when 
constructing the gravity dual of a confining theory, the KK scale and Aged are of the same 
order and cannot be separated within a weakly curved framework. 

The discussion above should be compared with what happens in unfiavored backgrounds, 
like (j2.H) - (j2.8p . the KS background jB] or the ones we present in section |H1 where there is no 
possibility of extremizing the function L{po). However, notice that the behavior found above 
is very similar to the one observed in the paper jHE] for a gravity dual of YM* without flavors. 
Such a background, like ours, presents a good singularity in the IR and the authors of [SB] 
considered the Wilson loop behavior unphysical and related to the fact that the string wanted 
actually to fall in the singularity. Nevertheless, there is an important difference between the 
two models since (and therefore the QCD string tension) never vanishes in our setup. 



5.3 Instanton Action 

Let us compute the action of an instanton in the field theory. We will propose that a gauge 
theory instanton should be thought of as a euclidean Dl-brane wrapping some contractible 
two-cycle in the geometry ()4.H1 . The action of the Dl-brane is: 

S = Ti f d^ze-^^/d^g + T,t f C(2) = ^ + tOsqcd (5.11) 

J J 9 sqcd 
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where we have compared with the action of the instanton in order to get an expression for 
the four-dimensional gauge couphng and theta angle. Since dF^^,) 7^ 0, it is not possible to 
define a two- form potential such that F(3) = dC{2). But if we consider the subspace given by 
6 = 6, ip = 27r — (p (and take also large constant p) we can write, from ()4.12|) : 

F(3)e// = ^(2 - x) sm6d6 A dip A d^j (5.12) 

and therefore: 

/V 

C(2)e// = (V^ - V'o) - ^) sm6d6 A dif (5.13) 



We now pick the two-cycle defined by 

6 = 6, (p = 2n-ip), ^ = {2n + l)n (5.14) 



So we can identify 
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6sqcd = ^(2iV,-iV;)(^-^o), , (5.15) 



Ssqcd 4 



— - e^^' + ^ia-lf (5.16) 



It is quite nice to see that the coupling constant satisfies the confining behavior we expect 
for Nf < 2Nc. Indeed, plugging in the expansions for the functions near the origin ()4.21|) . 
we see that even though the solution is singular, this singularity is not causing problems and 
the coupling grows at low energies. Conversely, in the UV, using the expansion in ()4.19|1 . 
we observe that the coupling is asymptotically small. We will further discuss the quantities 
written in ()5.15|) . ()5.16|) in sections EiH and EUl respectively. 



5.4 U{1)r Breaking 

In the previous subsection, we have identified a contractible two-cycle on which we wrapped 
a euclidean Dl-brane that we identified as an instanton. Since this cycle contracts to zero 
near the origin of the coordinates, the partition function for this euclidean brane should be 
equal to one, up to a phase. For the original presentation of this argument see jlU]. The 
partition function of this instantonic string is 




(5.17) 



^^We remind the reader that we have used that the tension of the Dl-branc is Ti = ^r-^, — and we are 

^ lira Qs 

using a' — gs = ^ along our computations. 
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where S2 is the cycle ()5.14|) and (...) represents the Born-Infeld part of the action, see ()5.11|) . 
Under a shift in the i?-symmetry ip —>■ + 2e (we write the factor of 2 since ip G [0, 4tt) in 
order to have e G [0, 2tt)) we will have 

Z Z[0] exp[te{2N^ - Nf)] (5.18) 

so, asking for invariance of Z imposes that (this is equivalent to requiring that the Osqcd 
defined in (I5.15p is unchanged modulo 27r under shifts of ip): 

^'^^ :^ = l,...,2Nc- Nf (5.19) 



2Nr-N 



f 

So we see that quantum effects select a discrete subset Z2Nc-Nj out of the possible rotations 
of the angle ip associated with the i?-symmetry, that is consequently broken according to 
U{1)r Z2Nc-Nf Z2, with the last step in the chain of breakings due to the fact that the 
full background is only invariant under ip ip + 27T. This last step in breaking is understood 
as the formation of a condensate, to which the functions a{p), b{p) should be dual to ^Tj . 

Notice that when Nf = 2Nc the breaking above does not happen and the U{1)r symmetry 
preservation indicates a theory with higher invariance, in this case, scale invariance: the 
U{1)fi anomaly is in the same anomaly multiplet as the beta function so, when Nf = 2Nc 
the coupling in ()5.1f)|l should not run. We study this point in more detail in section IHl 

5.5 Beta Function 

Let us start by summarizing the result of the wilsonian beta function for SQCD 

/5. = -|f^(3iVe-iV^(l-7o)), (5.20) 

where 70 is the anomalous dimension of the quark (and anti-quark) superfields, while for the 
coupling of the quartic superpotential we have 

/?;,~(l + 27o). (5.21) 

An approach to the calculation of the beta function of SYM from a gravity dual is presented 
in Without entering into the many possible discussions about the validity of this 

computation (more on this in section IHTSj) . let us just repeat the steps in jSH] for our case. 
The coupling is defined in ()5.1(j|l and we define the radius-energy relation near the UV of 
the field theory (at large values of p) as 

(A)3 _ a(p) ~ b(p) ~ e^20 (5 22) 
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where with the symbol ~ above we mean that all these definitions, when taken at large values 
of the radial coordinate, give the same result. The idea of the first two definitions is that 
these functions should be dual to the gaugino condensate or other condensate as proposed 
in IH]. On the other hand, one might consider that the exponential of the dilaton is related 
to the strong coupling scale as proposed in [3ll2j. Considering only the large p expansion 
of the coupling defined in ()5.16p and the functions in ()5.22|) , we have 

^~iV,(l-|)p + ...., logf~|p (5.23) 

dsqcd ^ A J 



So, a straightforward computation leads us to 
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a = = -^(2iV, - Nj) + ... (5.24) 

a log f 327r^ 

where the last equality should be understood as valid only in an expansion for large values 
of p. This shows that in the regime described by the large p region of the background, 
the anomalous dimension of the quark superfield is 70 = —1/2 (and therefore Ph = 0). 
Once again, this computation should be severely criticized regarding the approximations and 
assumptions that are done (see section |H1 for a different sort of criticism to this approach), 
but we just wanted to spell it here to show that the coefficient 2Nc — Nf appears. This 
also reinforces the fact that for the case in which Nf = 2Nc something special is happening 
(and we know the solution changes qualitatively for that case). It would be interesting to 
understand what happens when Nf > 2Nc- 



5.6 Domain Walls 

The domain walls of A/" = 1 SYM associated to the IR spontaneous breaking of the R- 
symmetry 7L2Nc ^TLi correspond to D5-branes wrapping the finite 5''^ at p = [7]. Their 
associated worldvolume ^-symmetry matrix is: 

x^^xxx^viz^* = e (5.25) 

As pointed out in jlH], this projection commutes with those in ()B.4|) if and only if ^ = 0, i.e, 
when the brane is placed at p = 0. So from the brane point of view, one finds, consistently, 
that these objects are one half BPS. 

In the fiavored theory we are dealing with, we have argued that there is a also a spontaneous 
IR breaking 21 2Nc-Nf — > ^ 2 and therefore one would also expect to have domain walls. It is 
easy to see, using the expansion ()4.21|) and the expressions in appendix |Bl that it is also true 
in our fiavored case that limp^o-^ = 0, so the same kind of embedding at p = preserves 
again half of the supersjTiimetry. Notice also that the action for the domain wall D5-brane 
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wrapping the three-cycle parameterized by 9, (p, ip at constant p is: 

En, oc T5 e"^^'^^' oc (^.iVj^iV, e^^° (1 + 2p2 + ...) (5.26) 

6Ci 

and it has indeed a minimum at p = (we have used eq 1)4.211) and restored the dependence on 
Qs). Also, despite e^^ and e^'^ tend to infinity and zero respectively, the relevant combination 
()5.26|) goes to a constant at p = 0, rendering a finite tension for the domain wall (its precise 
value would depend on the integration constants ci, C2 and could be used to fix them). Thus, 
we would like to stress that, even if at p = we are precisely on top of the curvature 
singularity, it is appealing to find sensible results from these simple computations. 



5.7 Seiberg Duality 

A remarkable feature of many Af = 1 gauge theories that include matter transforming in the 
fundamental representation of the gauge group (or even in the bifundamental representation 
for a quiver gauge theory) is Seiberg duality In its simplest form, this can be stated as 
saying that an A/" = 1 SU (Ac) gauge theory with Nf quarks and an A/" = 1 SU{Nf — N^) 
gauge theory with Nf quarks and a gauge-singlet meson have the same IR physics. When 
the theory contains a quartic superpotential for the fundamental fields, the duality is a bit 
more complicated |45j but nonetheless it retains the same relations between the number of 
colors and fiavors of the two Seiberg-dual theories as in the simpler case. 

In this section we will argue that Seiberg duality appears in our construction in a remark- 
ably simple and elegant way: it corresponds to switching between two alternative ways of 
writing the background^^ . 

The central point on which the realization of Seiberg duality in our model is based, is that 
the internal space of our background has S*^ x topology, but the choice to use either the 
space spanned by the pair {9, ip) or by {9, (p) as the 5*^ base of the 5*^ is (at least at first 
sight) arbitrary. In fact, let us take the metric and F(3) fiux to be written as, from ()4.9|) . 

ds^ = 62 ( rfxi 3 + dr^ + e^'\zj + eg) + — {{coi + a(r)ei)^ + (^2 + a{r)z2f) + —ul j 



(5.27) 
and from ()0)) and 

E{3) = ^ {-oJi A CJ2 A CJ3 - b{r) (ei A (1^2 - ^2 A cui) A LJ3+ 



-h'{r)dr A (ei A cji + e2 A CJ2) - (1 - x)ti A 62 A u^} 



^^A first hint of what we wiU develop for the non-abehan background can be easily gained in the simpler 
abelian case by looking at the way the value of the integral of F(3) (|3.18|l changes when it is evaluated over 
the 3-sphere {^p,9,ip) or (ip,6,ip). 
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where the functions 0, a, g, h and k are numerically evaluated in section 14 ■4[ b is given 
by 

and for convenience we have defined the following 1-forms 



(5.30) 



ei 


= de 


ei 


= de 




= — sin 6 dip 


h 


= — sin 6 dip 


UJi 


= cos ip ti — sin ip t2 




= cos ipii — sin ip C2 


UJ2 


= — smip ti — cos ip Z2 


UJ2 


= —sinipci — cos ip t2 


Us 


= dip + cos 6 dip 




= dip + cos ^ dip 


0)3 


= + cos 6 dip = + cos 







We know that the number of colors of the dual theory can be evaluated by integrating the 
three-form flux over the internal spanned by the angles 6, ip, ip, a.s in (j3.4|l 

NcT, = ^ / i^(3) (5.31) 
^^(10) J53 

whereas the number Nf of flavors in the gauge theory is determined by the failing of the 
Bianchi identity for F(3), equation ()4.1()j] . 

It is easy to show that the 1-forms we defined in ()5.30|) satisfy the following identities'^ 

UjI + UJ2 = tl + tl ul + ul = tl+ l\ 
tl Ui + e2 (1^2 = tJi tl + UJ2 t2 

cui A = — ^1 A e2 cJi A CJ2 = — ei A C2 (5.32) 

ei A c<j2 — ^2 A cji = — (ei A 0)2 — e2 A uJi) ti /\ uJi + t2 /\ uj2 = — (ei A cDi + e2 A CJ2) 

We can use these expressions to rewrite the metric and F^^,) iii a new form. It is very 
important to notice that even though these two quantities will look quite different from 
above, what we do is only a rewriting of ()5.27p and ()5.28|) . For the metric we obtain then 

ds'^ = 6 2 ( dxl^ + dr^ + e^^{t\ + e^) + — [{uji + d{r)tiY + (^2 + a(r)e2)^) + —CjI \ 

(5.33) 

and for the three-form flux 

-^(3) = ^ {-^1 A CJ2 A el's - h{r){pti AiJ2-h/\ t^i) A LJ3+ 

-h'{r)dr A (ei A cji + e2 A UJ2) - (1 - x)li A 62 A c^s} 



^■^We thank Agostino Butti for bringing these identities to our attention while describing to us the real- 
ization of the SU{2) X SU{2) isometry of the baryonic branch backgrounds of ,24) . 
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where in fl5.33|) we have defined 

/ ^2^2g 



e^' = — I 1 



4 V Ae^^' + a?e^\ 



a = a 



and in ()5.34|1 we have set 

N, = Nf-N, x=-^ bir) = P^ (5.36) 

X — 1 1 — X 

in order for the new expressions of the metric and the three-form fiux to have again the 
same apparent form as in the original versions (|5.27p and (j5.28|) . By using the definition of 
X = Nf/Nc and the expression for Nc, it is immediate to show that Nf = Nf. In the same 
way, it follows from ()5.36j) and ()5.29j) that 

6(p) = (2-x)-^ (5.37) 
smh2p 

that is, exactly the same form as in ()5.29p with the new value of the Nf/Nc ratio. 

Because of the way the metric looks now (j5.33|) . it seems natural to identify as the three- 
sphere of the internal manifold, the one spanned by 6, (p and ip. Therefore evaluating the 
integral of F(3) over this S*^, we find that the number of colors of the field theory dual to the 
background (|5.33p - (j5.34p is N^. = Nf — whereas the number of fiavors is still Nf. This is 
very interesting because this is exactly the way the rank of the gauge group changes under 
Seiberg duality. 

Let us now compare the geometries in the IR, where Seiberg duality tells us that the two 
field theories are the same. By using the small p expansions of section we can derive the 
behavior of h, g and a near the origin 

e2^ = e2'^--(2-x)p^ + ... 
3 

^2-g^^29^^^2-x)p' + ... (^-38) 

a = a — 2cip^ + . . . 

and remarkably we see that the newly defined functions correspond to their original homonyms 
up to the third term in the small p series expansions. This is the holographic dual of the 
statement that the Seiberg pair theories have the same IR behavior. 

This result might look a little puzzling at first sight: in ()5.33|) and ()5.34|) we have only 
written in a different way the metric and three-form fiux but at the same time claim that 
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the background written in this way corresponds to a different theory from the one that is 
dual to the (same) original solution ()5.27|1 and ()5.28j) . We believe the right way to look at 
this puzzle is in terms of the dictionary we use to translate the geometry data into field 
theory predictions. So in a sense, we have two different dictionaries we can use on the same 
background, corresponding to the two ways of identifying the inside the internal space. 

This is a remarkable result: just by rewriting the flavored background in a different form, 
we can capture the most relevant features of Seiberg duality for theories with flavors. In 
fact in a simple and elegant way we can reproduce the change in the rank of the color gauge 
group while leaving the number of flavors untouched, and we can show that the IR dynamics 
of the two theories is exactly the same thanks to ()5.38j) . It would be very interesting to 
study how information about the Seiberg-dual field theory is encoded into string theory in 
this setup, and which predictions this allows us to make. 



6 Features of the Nf = 2Nc Solution 

We now discuss the gauge theory interpretation of the particularly simple Nf = 2Nc solution 
{or X = 2 in the notation introduced above), presented in section see eq. fl4.23j) . For 
this solution (and only for this one), the effective four-dimensional Yang-Mills coupling, see 
()5.16|) . is constant (notice that the couphng is constant regardless of the cycle used to define 

It is also the only solution for which the t/(l)_R symmetry is preserved. There is not anoma- 
lous UV breaking (we can apply the discussion of section EH eq (j5.18j) ). while there is not 
spontaneous IR breaking since a = b = 0. We would therefore not expect to have domain 
walls and this is in fact the case since its tension as defined in ()5.26|) would vanish (e^"^ 
when r — > — oo). We think that these two features are signaling an IR conformal fixed point 
in the effective four-dimensional theory. 

That this happens only for Nf = 2Nc is in precise agreement with the field theory in- 
terpretation of section 15.11 Notice that, unlike in pure A/" = 1 SQCD, there should not 
be a conformal window for Nf < 2Nc. In fact, having an IR effective quartic superpo- 
tential for the quarks fixes their R-charge to |. The NSVZ beta-function, which can 
be written in terms of the (IR) R-charges of the fields can consequently only vanish when 
Nf = 2Nc. This can also be seen using the wilsonian beta-functions for the gauge coupling 
Pg ~ {3Nc — Nf^l — 'jo)) and the quartic coupling ~ 1 + 270 which vanish simultaneously if 
and only if 70 = — | and Nf = 2Nc (70 is the anomalous dimension of the quark superfield). 

Nevertheless, we cannot claim that the solution (j4.23|) is a faithful dual to a conformal 
four-dimensional theory. Since there is no AdS^ space, the conformal group 5*0(2,4) does 
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not show up in the gravity side and also there is no superconformal enhancement of the 
supersymmetry. We beheve that the explanation resides on the fact that we start with 
a six-dimensional theory living on the D5-branes compactified with a twist on an 5*^ and 
then add flavor branes. The four-dimensional theory never completely decouples from the 
Kaluza-Klein degrees of freedom and this spoils the would-be conformal symmetry. Indeed, 
even though the parameters and couplings might be tuned so that the massless modes are 
in a "conformal phase", the KK modes spoil this symmetry and the background reflects it 
by not allowing an AdS^ to appear. To add to this interpretation we can, for instance, make 
a dilatation of the four dimensional spacetime Xi Xxi, that can be compensated (we take 
string frame) by a shift r — r — 2 log A, which should be regarded as a rescaling of the energy 
scale (see ()5.22|l ). This is not an isometry of the metric since the internal space gets rescaled, 
thus rescaling the KK scale. However, we find appealing that the geometry is not deformed 
in any other way, since the functions h, g, k, a do not depend on r. Notice that this happens 
only for this simple Nf = 2Nc solution. 

In any case, one could not expect to find a weakly curved AdS^ since this would imply that 
the coefficients of the conformal anomaly are equal a = c at leading order in Nc [IE]- Since 
our theory has Nf ~ one has a 7^ c at leading order 07], ruling out the possibility of 
having a faithful supergravity dual. Nevertheless, we think it is very nice that some features 
can be matched in this setup. 

Finally, let us conjecture on the gauge theory meaning of the parameter ^. A/" = 1 
SQCD with Nf = 2Nc and a quartic superpotential has a one complex-dimensional family 
of marginal deformations in the two complex-dimensional parameter space of the gauge 
coupling and of the coupling of the quartic term (see The parameter ^ changes the 

volume of the spheres ()4.22|) so it is natural to think that it would change the masses of the 
adjoints coming from the KK modes and, indeed, modify the value of the coupling at the 
"conformal fixed point" ()fi.H9j) . Thus, it may be a parameter signaling this conformal line. 

Let us see this more precisely. In this geometry, we could define couplings (not necessarily 
the SQCD coupling in ()(j.39|l ). to be proportional to the inverse volume of two-cycles 

1 Vol S2 



94^ 9I 



(6.40) 



Here = a'ggN is the coupling of the little string theory. So, different definitions would 
imply the choice of different two-cycles in the geometry, like 

1 _ Vois\e,^) 1 _ Vois\l^) 1 _ Vois\e,l^,^) 

^ 9i 9^ 9i 9;qcd 9i 

Let us be more precise and define the cycles above. The first two cycles are given by the 
effective geometry obtained from ()4.23|) when imposing that all the coordinates are constant 
except {0,(p), or {6,ip) in the case of S'^{6,(f) and S'^{6,Lp) respectively. The last cycle is 
defined like in ()5.14j) and gives the coupling in ()(j.39|l . 
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So, computing these "couplings" one gets 



^Eh, ^ ' 



-En, 



1 4 

-0 r— ^ (6.42) 

Where i?[7r, /c] is the second kind Elhptic function. Now, let us observe something quite 
nice. If we interchange the cycles S^{9,ip) ^ S^{9,ip), this is equivalent to the interchange 

^ (4 — ,^), or ^ ^9- The Seiberg-like duality defined in section IKTI does precisely this 
job. So, we have a line of fixed points, parametrized by C, (solutions where the two-spheres 
have constant radius), which is reflected under Seiberg duality. This is what is known to 
happen in a field theory like the one discussed in section that is A/" = 1 SQCD plus a 
quartic potential, when the theory is at its conformal point {Nf = 2Nc, and the anomalous 
dimension of the quark superfields is 70 = — |). For a review of these results, see Strassler's 
beautiful lecture notes |45]. Indeed, following jlSj, we can see that the line of fixed points 
has a self-dual point (under Seiberg duality), that in our case is ^ = 2. Then the duality, in 
our case the interchange ^ (4 — ,^) or the interchange of two cycles, interchanges points 
on the fixed line that lie on different sides of the self-dual point. We can see that this is 
happening in eqs. ()6.42p . for a holomorphic coupling that should be a combination of gsqcd 
and the quartic coupling. 

This suggests that as for the coupling ()6.39|) . it might be possible to find another two-cycle 
whose volume is inversely proportional to the quartic coupling. This is a very nice check 
that points to the consistency of our approach and interpretation. 

Finally, we would like to stress two very interesting facts. First, let us compute the Wilson 
loop (see section eq. ()5.10p ) for this particular case. It can be seen that for any value 
of one obtains that the quark antiquark pair has energy Egg = and that the separation 
of the pair is always L = iry/NcQsO:'. Following the interpretation of section 15. 2[ below 
this length this computation is unphysical due to the UV completion. For L > Hy/Ncgsd' , 
the configuration in which the qq are joined by a string is never preferred to having two 
independent particles (notice that in this case, the QCD string tension vanishes gtt{r = 
—00) = 0). We interpret this observation as the existence of total screening, the "mesons" 
have zero energy, in agreement with the fact that this theory does not dynamically generate 
a scale. 

The second observation concerns the black hole solution in ()4.26|) . We can compute the 
viscosity of the dual quark-gluon plasma following the approach developed in jlHj. Since our 
background ()4.26|) satisfies the hypothesis of the theorem in jlH] , we can see that the relation 
between the shear viscosity and the temperature is of the predicted form V = It may 
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be interesting to study other hydrodynamical quantities using the dual solution ()4.26|) . This 
is a very simple background dual to a field theory with adjoints, fundamentals and at finite 
temperature. 

Finally, to close this section, we would like to call the reader's attention to the solution 
presented in appendix O There we described a numerical solution for the case Nf = 2Nc, 
that appears after a careful analysis of the BPS eqs. Interesting features of this solution are 
that it has a behavior similar to the fiavored case (with Nf ^ 2A^c), asymptotes in the UV to 
the simple solution presented in eq. ()4.23p and the fact that the solution is Seiberg-duality 
invariant. We do not know how to interpret this new solution from a dual field theory 
viewpoint. It seems to behave as a fiow from a fixed point where a relevant operator is being 
inserted, breaks "conformality" and goes down to an IR theory similar to the one described 
by the solutions with N f < 2Nc. We leave this interesting problem for future study. 

7 An Alternative Approach: Flavors and Fluxes 

As we have seen in section 14.21 the addition of many fiavor branes to the MN background 
is encoded in the non-closure of the three- form RR fiux F(3). When the fiavor branes are 
smeared over their transverse directions, the failing of the Bianchi identity is proportional 
to the volume form of the transverse S"^ x S"^ 

Nf ~ 

dF3 = ^ sine sine de A d(f A de A d(f. (7.1) 

In sections El and E] we shown that this effect gives rise to a background which has a uniform 
distribution of F(3) charge along the two two-spheres. As we have shown in the preceding 
sections this procedure works very nicely and gives rise to a remarkable matching between 
gauge theory expectations and string predictions. 

Nonetheless, we think it could also be very interesting to approach the same problem 
from a complementary approach. In this section we sketch the general idea and some basic 
features of this possible new scenario, but setting up and solving the whole problem would 
require a much deeper development and thorough discussion, which we leave as a possibility 
for the future. 

In this new approach, the failing of the Bianchi identity for the three-form F(3) is achieved 
through non-trivial fiuxes only. Indeed, we know from supergravity that the presence of an 
NS three-form background deforms the definition of the field strengths of the RR fiuxes. Let 
us suppose to turn on some axion C(o) and NS fiux H = dB(^2), then the field strength F(3) 
is defined by 

Fs = dC^2) - C^o)H (7.2) 
and we see that if the axion is non-constant we can have a non-trivial Bianchi identity for 
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^(3) 

dFs = -dC^o) = A H (7.3) 

Notice that in this case the decomposition of -F(3) in an exact and non-exact part is natural, 
and therefore the definition of C(2) is unambiguous (up to gauge transformations, obviously). 

In the preceding sections we have always assumed that the susy transformations of the 
flavored background are exactly the same as those of the MN solution [3 . This is due 
to the fact that we are adding to the MN solution an (infinite) set of BPS branes, which 
by definition preserve the same supersymmetries as the original background. But when 
we look at the flavoring of the background from the flux side, this could fail, and in fact, 
this is what one usually expects when turning on new fluxes in a background. The only 
condition we think it is still sensible to impose, is that the background can accommodate 
supersymmetrically embedded branes. 

Regarding the geometrical properties of the six-dimensional internal manifold. There are 
two possibilities to build backgrounds with fluxes that preserve four-dimensional A/" = 1 
supersymmetry: SU (3)-structure manifolds [20] and a particular class of S'f/(2)-structure 
manifolds jST]. The unflavored MN background the internal manifold has S'f/(3)-structure 
inO] and since we can most likely think of the flavored background as a deformation of the 
MN solution, it seems reasonable to assume that its internal part is an S'[/(3)-structure 
manifold as well. 

Starting from these two assumptions, we can write the spinor ansatz for the supersymmetry 
transformations. As usual for solutions with four- dimensional Lorentz invariance, we split 
the IIB ten-dimensional spinors in their four- dimensional and internal six-dimensional parts 



£2 = C+ ® ^+ + C- ® ^- 

where r/^ and rf are spinors of the internal manifold. Notice that r/!^ = (?]+)*. Since S'f/(3)- 
structures allow only for one globally-defined spinor r/, we need to take r)^ and rf to be 
proportional 

ri\ = ar]+ rjl = br]+ (7.5) 

A necessary condition to have consistent supersymmetric embeddings of branes in an 
S'f/(3)-structure manifold, is that the two spinors r]^ and rj'^ have the same norm, = 
II^^P |S2]; which requires that the functions a and b in ()7.5p satisfy 

^ = e'^ ^ = e'^ (7.6) 

with and r real functions of the internal coordinates. For the smeared configuration we 
consider in this paper, symmetry imposes that they only depend on the radius of the internal 



41 



manifold. Notice that the spinor t] is defined up to a phase, and therefore r in ()7.6|) can be 
arbitrarily gauge fixed by this symmetry 50^- On the contrary (p does not depend on the 
gauge choice. 

Instead of a and b we can define a and jS 

a = a + ib /3 = a — ib (7.7) 

in such a way that when we define e = ei + ze2, we get 

e = aC+0V+ + PC-®V- 

In terms of a and /3 the first condition in ()7.6|) reads 

^ = -^cot| (7.9) 

At this point, one should write a smart ansatz for the metric and background fiuxes, and 
then the SU (3)-structure techniques developed in [SUl HH] would allow to derive first order 
differential equations describing the solution. This is beyond the scope of the present work. 



8 The Unflavored Case Nj = 0: Deformed Solutions 

This section might be read almost independently of the rest of the paper and addresses the 
problem of finding gravity duals to different UV completions to A/" = 1 SYM from the one 
described in section 2. Concretely, we will deal with the setup presented in section 14.11 
These solutions have appeared, in the context of this work, because in the procedure we 
devised to add flavors, we flrst deformed the original (unflavored) background ()2.H) - ()2.8|) . 
see the details in section |3] As explained in section these backgrounds are an interesting 
subcase of the ansatz studied in pij . 

Below, we will study a little more the interesting solutions of the unflavored system (j4.H) - 
fl4.8|) . In section IHUl we describe the backgrounds arising from the study of these equations 
and in sections I8.2H8.8I we discuss their gauge theory interpretation. 



8.1 Description of the Solutions 

Equations ()4.6|1 have a one-parameter family of regular solutions. In fact, for small values 
of the radial coordinate, we flnd a very well controlled expansion in terms of a parameter /z 
taking values in the interval (—2,-2/3]: 

2 , 2i / 4 20 + 36,^ + V 2 



a(rt = l+,p^ + .... e- = iV.(— -^5p^p^ + ..„) . (8,1) 
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By inserting it in ()4.7|) . ()4.8|) we find, for tlie otlier functions: 



In the case in which the parameter /x = —2/3 we reobtain the known solution ()2.H) - ()2.8|) . 
that is: 

e^^ = Nc (8.3) 



sinh 2p 

and in the case fi = —2 (and = 0) we get four- dimensional Minkowski times the deformed 
conifold (see appendix [nj. For p — > oo (except in the particular /i = — | case), the solutions 
asymptote to the deformed conifold metric (for details, see appendix [HI), i-e very far away 
from the branes their effect becomes negligible and the background asymptotes to a Ricci 
flat geometry. 

In the figure ini we present numerical plots of several of the functions defining the geometry 
for different values of /x. Notice in particular that, unlike the usual solution ()2.H) - ()2.8|) . here 
the dilaton does not diverge^^. 



8.2 Gauge Theory Analysis 

All the deformed solutions presented in section 18.11 are non-singular duals to field theories 
with field content being a non-abelian gauge field, a Majorana fermion and a set of massive 
adjoint scalar multiplets. These field theories in principle confine since the warp factor in 
the metric 1)4.11) goes to a constant value near p = 0. This will be analyzed in more detail in 
section ISTI Regarding the string theory picture, this is a compactification of Nc D5-branes 
on a two-sphere, that preserves Af = 1 SUSY. As studied in many places, and recently 
reviewed in detail in (20], this leads to a field theory whose lagrangian is the one of A/" = 1 
SYM coupled to a set of massive scalar multiplets. We will assume this picture below, and 
will argue that the dynamics of the massive modes is related to the parameter that labels 
the family of solutions. 

Contrary to what one might think, this is not the "brane side" of the solution ()2.H) (for 
the case of D6-branes this "brane side" solution was found in jMI)- If these new solutions 
were the "brane-side" , one should see a resolved conifold and the metric structure that 
characterizes the D5-branes, near the origin of the space. One might want to interpret this 
new solution as the "non-near-horizon" version of the solution of D5-branes wrapping a two- 
cycle inside the conifold. This must be at least partially correct, since at p — oo one goes to 
Ricci flat geometry and the effect of the fluxes disappears (notice also that a metric like 1)4.1)) 

^''A usual criticism to solutions like (|2.1(l - (|2.8)l is that since the dilaton diverges for large values of the 
radial coordinate, one must S-dualize and we should be dealing with NS5 branes, hence, the UV completion 
is a little string theory, which is an unconventional completion for a field theory. 
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Figure 3: We plot the different functions for different values of the parameter /x. The thick 
solid line corresponds to /i = — | (the usual MN case), the solid thin lines correspond, 
respectively, to /i = —.68, —1, —1.5 and the dashed hne is a deformed conifold. 

cannot be obtained by uplifting from 7d gauged sugra in the same way as ()2.1|) ). However, 
it cannot be the whole story since, as we will see, the IR physics is also modified within the 
family of solutions. 

The proposal we want to put forward is that different members of this family (that can be 
thought of as different values of fi in the window (—2, —2/3]) differ in the superpotential that 
affects the dynamics and masses of the KK modes. At the same time, it is also possible that 
different members of this family differ in VEV's for some operator. Then, the UV completion 
of the theory changes, modifying the large p behavior of the background (this is in analogy 
with turning on an irrelevant operator in AdS^ x S^). But this change in the UV dynamics 
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also alters the infrared (in analogy with what would be a dangerously irrelevant operator). 

Indeed, we know that all these infinite solutions preserve A/" = 1 SUSY by virtue of a 
twisting procedure in the gauge theory. It should be quite interesting to analyze how the 
developments in [201 apply to this new set of configurations. It is also interesting to notice 
that any expansion that starts out of the window (—2, —2/3] for the parameter /i mentioned 
above generates pathological solutions. This suggests that the dynamics of the KK modes is 
such that it shows some instability or similar sick behavior for solutions that are outside this 
window (as the supergravity solution becomes singular) Besides, we stress that the values 
fi = — 2,yU = —2/3 correspond to conifolds with maximal and zero deformation parameter. 
So, from the gravity side the existence of this window is clear, see appendix O 

A first check of our proposal comes from constructing new backgrounds doing solution- 
generating transformations on these //-family of solutions, that have a ?7(1) x f/(l) isometry. 
These isometrics are global symmetries of the KK part of the spectrum, hence, constructing a 
new background by SL{3, R) transformations based on those t/(l)'s (following, for example, 
the techniques developed in jHSI), will only affect the dynamics of the KK modes. This 
idea was proposed in jTHl HO] to be the way to distinguish between a pure SYM effect 
and the dynamics of the KK modes influencing some observable. When we perform the 
transformation on this family of backgrounds (we will not report here the formulas, since 
they are quite lengthy), we observe a mixing between the parameter of deformation 7 and 
the parameter /i (changing the dynamics of the KK modes). 

In the following, we will study a sample of gauge theory observables that range from the 
possible existence of massless axionic glueballs, k-string tensions, annulons and PP-waves 
spectrum, Wilson loops and beta function. Technical details can be found in appendix El 
These examples will support the validity of the proposal that we have made above. 

8.3 Axionic Massless Glueballs 

In the paper 56J, the authors studied the possibility of flnding a massless glueball in the 

exact solution ()2.1|) - ()2.8p . They found that this glueball was not normalizable, which led 

them not to consider it to be a physical excitation. This is in contrast to what happens for 

the Klebanov-Strassler solution or a non-supersymmetric deformation of it in which there 

are, in fact, physical massless glueballs [3^1 137j. 

In the exact ()2.Hl - ()2.8j) solution {fi = — |), the dilaton diverges at large values of the radial 

coordinate and this is the root of the non-normalizability of this mode. In this new set of 

solutions, we have that for /i 7^ —2/3 the dilaton is bounded at inflnity, see flgureEl It is 

natural, therefore, to wonder whether this family of solutions might have a massless glueball. 

Let us analyze this in some detail. In principle, we do not expect this massless glueball to 

^^One might speculate with the existence of a potential term in the KK lagrangian of the form ($ denotes 
a generic KK mode) V = - (a* + 2)(m + 2/3)$2" + ... that captures this behavior. 
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exist, since there is not in this model any apparent U{1)b symmetry that is spontaneously 
broken, leading to a massless supermultiplet of composites. 

In the paper the fluctuated eqs of motion were studied in detail. Let us assume that 
the background fields vary according to 

9pu Qpu + (^h^y, (j) ^ (j) + e6(f), F^yp F^yp + eSF^^p (8.4) 

So, keeping only linear order in the parameter e, we get eqs for the fluctuated fields that we 
report in detail in appendix |El One can easily see that a solution to these eqs is the one 
found by the authors of 

5F3 = *4rfa, V = = 0; = (8.5) 

Once again, the point in the paper j56j, is that this fluctuation, due to the diverging dilaton 
is not normalizable, so they reject the excitation as a physical state. We might wonder what 
will happen now that we have a new family of solutions with a bounded dilaton. The norm 
of the fluctuation is 

poo 

In this case the dilaton does not diverge, but the functions g, h, k are unbounded making again 
the mode under consideration not normalizable. Even in this family of solutions, therefore, 
we do not seem to have this massless glueball. There might be another combination of fields 
that could play that role, but this does not seem likely, due to the fact that this family 
of theories (if they differ from ()2.1|) - ()2.8|) in the masses and dynamics for the KK- modes 
and possibly condensates) should not have a U{1) symmetry that is spontaneously broken 
leading to a massless excitation. 

8.4 Confining /c-Strings 

One very interesting problem in SU (A^)-QCD and similar theories is the study of flux tubes 
induced by color sources in higher representations. If the source has k fundamental indexes, 
the flux tube is called the /c-string. The law for the tensions of these /c-strings has been 
subject of interest, from the lattice viewpoint and also from a more formal side. 

In the very nice paper Herzog and Klebanov tackled this problem using AdS/CFT-like 
dualities for confining theories. They considered a confining string (connecting quarks in the 
antisymmetric representation) as a D3-brane wrapping a cycle near the IR of the geometry. 
They found the tensions of these strings. Let us briefly summarize their approach and apply 
it to our cases of interest. This can be done easily with the expansion written in ()8.H) - ()8.2|1 . 
since we have to look at the geometry near p = 0. It is easy to see that our geometry (after 
the choice of the cycle 9 = 9,ip = 2t^ — (p^ip and a rescaling ip —>■ ip/2) turns exactly into the 
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geometry in eq.(12) of jSH]. So, near p = the metric induced on these D3-branes that wrap 
the cycle above and extend in [t, x) forming a confining string, reads 



dxl^ + e^'^^o) {d^^ + cos^ ^{de^ + sin^ ed^^)) . (8.7) 



Following the approach of Herzog and Klebanov, we can compute the tension of the q- 
confining string as the tension of a D3-brane wrapping the finite three-cycle above 



T 



The coefficient is 6 = ^-^p- = q^^- Only for /z = —2/3, that is for the exact solution ()8.3|) . 
we have 6 = 1. It is precisely this coefficient the one that propagates into formulas (12)-(18) 
of jnn]- Then, when minimizing the action of this confining string, our formulas will be 
exactly like the ones in . The tension of the confining string is minimal when 



r 2 • 4 , , / , sin2V' 7rg 2 
sm -w + (w 



1/2 



7/- = ^y^sin?/^ + ^, Tg = bsm^lj^/lTW^-Wco^ (8.9) 



The law that is conjectured to work for = 1 SYM (see for example jHO]) is 

Tfc _ sin(7rA;/A^c) 
~ sm{nk'/N,) 



(8.10) 



whereas the string tension in this family of models is not obeying a sine law, and only for 
the exact solution (// = —2/3) this behavior is reproduced. It is interesting to notice that 
observables like this seem to diverge for /i = —2, but we know that this point is out of our 
family. 

We mentioned above that we could interpret the parameter fi as different masses, and 
most likely also different parameters in the superpotential, for the KK modes. This matches 
nicely with the findings of the paper jHIj. The authors there have pointed out that if one 
considers N" = 1 SYM as broken M = 2 SYM, the sine law is not universal. The set of 
models we have resembles quite much these A/" = 2 — >■ A/" = 1 models (it is a little more 
subtle since we have a twisting). In the analysis done in jJU] for the Klebanov-Strassler (KS) 
model, it was found that confining strings there do not present a sine law behavior. In this 
aspect, the family of deformed solutions we are studying, behaves hke the KS solution. 



8.5 Rotating Strings 

We will study rotating strings in the same spirit as advocated by Gubser, Klebanov and 
Polyakov in the paper jHS]- These rotating strings should be dual to large-operators in the 
gauge theory (composed of a large number of adjoint fields). We will see that these strings 
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display characteristic energy-angular momentum relations and sample (for short strings) the 
different dynamics of the KK modes, that is, they depend on the particular value of fi. There 
are a number of solutions already found in [^3 that show this behavior. So, consider for 
example a string that rotates on S"^ and is stretched in the radial direction parametrized by 

t = KT, ^ = UT, 9 = -, p = p{a). (8.11) 

After writing the Nambu-Goto action, the energy and angular momentum can be computed 
(for a string that stretches from the origin to a position po) one gets jHSHHII 

PO g2(^+2fe 

^/^2^+2fc(^^^^^^^2^^2r^lp^ 

J~ / ^ ^ ' = (8.12) 

By doing an expansion near po ~ 0, that is by considering short strings, it is easy to see that 
the dispersion relations for this string depends on the parameter /i. This means that the 
operator dual to this string configuration contains a gauge invariant combination of fields, 
including the KK fields discussed above. In other words, these operators will have different 
energy-angular momentum relations depending on the member of the family of solutions we 
consider. 

Now, we would like to study string configurations whose energy and angular momentum 
expressions do not depend on the p parameter. The latter are more interesting than the 
solutions discussed above, since they should be dual to pure gauge theory operators made 
out of A^, A and derivatives; so we will be more detailed. For this, we consider a configuration 
for a bosonic string that rotates on "the gauge theory space" and on an internal direction 
(chosen so that it will give the string some i?-charge) and also stretches along the radial 
direction. We consider the ansatz (see appendix El for details) 

t = KT , X = R{(t) cos(ti;ir) , y = R{<j) sin(ci;ir) , 

p = p{a) , il) = 2uj2T (8.13) 

Let us take a look at some particularly interesting solutions. 




Rotating folded closed string 

Let us consider constant p. Then, one has to take p = uj2 = Q (notice that, at p = 0, we 
have dpk = dp(f> = 0) and: 

R = —cosuia (8.14) 
^^See also [HI] for the discussion of a similar configuration in the backgrounds of [7II?>K]. 
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This is dual to a high spin glueball. Notice that the slope of the Regge trajectory will depend 
on the tension of the string at the bottom of the geometry and therefore on e*^" . 

Pointlike strings rotating in 

Consider: 

R = const p = const (8.15) 

Then, we need uji = and: 

dpk = ^ p = 

e^^ool = 1 (8.16) 

These strings are dual to (long) R-charged operators . 

It is interesting to notice that when computing observables for these strings, like energy 
and angular momentum (i?-charges), they do not depend on the parameter p hence, they 
are the same for different members of the family. This can be interpreted as an indication^'' 
that these strings are a purely gauge theory effect and are dual to (large) operators that 
are gauge invariant combinations of A^, A. On the contrary, moving strings like the ones 
considered in ()8.11|1 - ()8.12|) (originally found in [ESI EI]), show dependence on the p parameter 
when computed in these new backgrounds. This is because they are dual to large operators 
containing KK modes, as established in jHl] (following the proposal of [IHl)- It would be 
interesting to recompute our solutions above in the deformed background proposed in [TH] . 



8.6 Penrose Limits and PP- waves 

This subsection is constructed on the work done in the very good paper . Let us summa- 
rize their findings. They studied a Penrose limit of the background ()2.H) - ()2.8|) on a geodesic 
that is mostly localized near p = 0. The limit gives a PP-wave geometry on which the string 
theory can be quantized. The authors of [Hni associated this geodesic on the string side with 
an operator on the gravity side that looks like a large chain composed of KK modes, that 
they called annulon. The spectrum of excitations of the annulon was put in correspondence 
with the spectrum of the string on the PP-wave geometry . We will follow their approach 
and take the same Penrose limit on our one-parameter family of solutions. The details of the 
computation are spelled out in appendix |El After following the approach in [H^ adapted to 
the present case, we get a plane wave background that looks 

g-2/c(o)~0(o)/2^2^^2 ^ _2dx+dx- + dxl + + dz^ + \dvl - (^ + ^yl)d(fi, 

F3 = {2dvi A dv2 + ^dzi A rfza) A da;+. (8.17) 

^^Notice that most probably //-independence is a necessary but not sufficient condition to ensure that the 
operator does not contain KK modes. 
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So, we see that we have four massless directions {x^, z) and four massive ones; V2 with mass^^ 
= 1 and y2 with = The excitations of this PP-wave are dual to operators 

composed out of KK modes and one can clearly see that the parameter yU is related to 
the mass of the KK modes. As we mentioned before, it may be the case that is also 
influencing the dynamics of the KK's or condensates, but this is erased in the Penrose limit 
and we only see the effect in the masses. So, following the treatment in jnO], we can compute 
the spectrum for strings on this family of pp-wave backgrounds and put the dynamics of this 
string in correspondence with the dynamics of the family of annulons. 



8.7 Wilson Loop 

We would like to analyze here a subtle point related to Wilson loops. We will follow the 
treatment developed in [SH] and carefully revised in jHE] • We will consider a probe brane that 
is wrapping the same 2-cycle as the branes that originate the geometry, but this probe brane 
is very far away from the origin (this probe was proven to be SUSY preserving in pH])- This 
leads to very massive fields that we associate with quarks, and we compute the expectation 
value of the Wilson loop as 

<W >r^ e^^« (8.18) 

where Sng is the Nambu-Goto action for a string in the background (jHH). After choosing a 
configuration for a string that extends in the coordinates x, p, t and is parametrized by the 
worldsheet coordinates ex, r according to x = ex, p = p{(j),t = r, the NG action will read 

SNG = ^j da^\g,^gu\ + \gugjp'^ = ^ / ^^/^"VT^'^VK (8-19) 

Notice that this action implies the existence of a conserved quantity, that we associate with 
F{po), where po is the point of maximal proximity of the probe string to the origin of space 
p = 0. The separation of the two very heavy quarks is given by L = 2F(po) fp^ ■j^^j==^^^dp, 

where pi is the distance where we put the probe brane (in principle one should take pi ^ 00 
to make the quarks very massive). 

So, following the treatment reviewed in [36|, one gets that the energy of the quark anti- 
quark pair is given by (normalized after we subtract the infinite mass of the non-dynamical 
quarks) 

E,^ = F{po)L + 2 ^^F^-F{poydp - 2 Gdp (8.20) 

Jpo ^ Jo 

The authors of proved a very nice theorem stating under which conditions the energy 
of the quark anti-quark leads to confinement. Some hypotheses are done on the behavior of 
the functions F{p) = e'^ and G(p) = e'^^'^ 

^^Notice that the mass of the two modes parametrized by V2 is not affected by the /i-deformation and is 
common a feature to all members of the family. 
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1) F{p) is analytic and near p = 0, F = F{0) + UkP^ 

2) G{p) is analytic and near p = 0, G = J2j ^jP^ 

3) F{p),G{p) are positive 

4) F'{p) is positive 

5) The integral J°° y^dp converges. 

We can see that the functions for our family of backgrounds behave, for small values of 
the radial coordinate, as 

So, one can check that hypothesis (l)-(4) above are satisfied but, unless p = —2/3, hypothesis 
(5) fails to be satisfied by the members of this family. This is quite remarkable, since it 
separates the exact solution occurring when p = —2/3 from the other members of the family 
(just like what happens when computing string tensions in subsection 18. 4j) . 

We want to argue nevertheless that all this family is dual to confining theories, since the 
SYM string tension will behave like T = F{0), but the correction to the leading behavior 
might be different from the one predicted in 13 6 j . It would be interesting to evaluate the 
integrals in ()8.2()|1 . to obtain the correction to the main confining behavior. 



8.8 Beta Function 

There is a computation of the N" = 1 SYM beta function done in In spite of being 

quite controversial (because some of the assumptions are not clear to be correct and one 
might criticize the regime of validity in which the computation is done), it works in a quite 
remarkable way, not only because the final result is correct, but also because it keeps working 
even when deforming the background in quite dramatic ways (like a non-commutative or 
dipole deformation [HZl CHI)- One way in which we might view this is that the authors of 
[39] have found a two-cycle (needed to define the coupling) that has very interesting and 
robust geometrical properties. What we will attempt in this subsection is to do again this 
calculation of the beta function. We will see that the result seems to change drastically from 
the one in (HH], but we will also point out an interesting subtlety regarding this point. 

Let us start by taking the expansion for large values of p given in eq. ()C.4|) . ()C.1|) . We 
see that the coupling constant, that is defined as the inverse volume of the cycle 6 = 6, 
if = 271 — (f, ip = 71 , is given by 

^ = ' (8.22) 
Using the fact that the gaugino condensate is related to the function b{p) [H] allows us to 
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write log(-E'/A) = 



— 1/3 log(6(p)). Then computing the beta function gives 









d{\og{fi/A)) \ 2 ) 




(c/log(/i/A)) 



-9 



(8.23) 



We see that this is not a running that we recognize as the one of a SYM theory. So, since 
this computation can be done with a member of the family very close to /i = —2/3 (for 
which the computation above reproduces the NSVZ result), we see that separating slightly 
from the exact solution ()8.3|) changes drastically the result. This may cast doubts on this 
computation being relevant for the SYM beta function. 

Nevertheless, there is an interesting fact that should not be overlooked. In Figured one 
can see that a coupling defined as in eq ()8.22|1 will behave, near the UV of the field theory, 
like 

^ (8.24) 



gm 



if we take the UV of the SYM field theory to be around the region bounded from above by 
the cut-off Auv, after which the stringy completion to the field theory should set in. With 
this caveat, following the procedure described in one gets the typical NSVZ result. This 
addresses one of the criticism that were raised to the papers j^^. We believe that this puts 
the whole result in a better context. 



0.2r 



0.1 





Figure 4: We plot the inverse 't Hooft coupling and e'^ for fi = —.67 (thin line) compared 
with the usual (thick line). Starting from some value of p (Auv), the UV 

completion sets in and the geometry asymptotes to a (Ricci flat) conifold. Auv decreases as 
/i decreases. 



9 Conclusions and Future Directions 

In this paper we have studied two very interesting problems. The first was to find supergrav- 
ity (string) backgrounds dual to minimally SUSY field theories in four dimensions containing 
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fundamental matter with Nf ~ Nf.. After writing the BPS eqs that characterize the back- 
grounds, we solved them asymptotically and showed numerically the smooth behavior of 
these solutions. Then, we studied numerous gauge theory aspects of these solutions, check- 
ing known results and making interesting predictions. We believe the material presented in 
this part of the paper may be quite useful in future approaches to non-perturbative SQCD- 
like theories. 

In a second part of the paper, we discussed a family of solutions dual to SYM plus some UV 
completion, that generalize the background of [7j. We studied the gauge theory aspects of 
this family of UV completions by analyzing the supergravity solutions. It is quite interesting 
to uncover the differences between gauge theory observables of different members of the 
family of solutions. 

We would like to comment now on possible future directions that this paper opens. Indeed, 
the first thing that comes to mind is that the formalism developed in the first part of the 
paper could be immediately generalized to find string duals to field theories with flavor in 
different number of dimensions and with different number of SUSY's, using backgrounds like 
the ones in |68j . 

It would be interesting, following for example the technology developed in [0^1 lEH]; to 
find non-supersymmetric deformations of the solutions we have presented here. This might 
prove a way to get information on QCD-like theories with an arbitrary number of colors and 
flavors. 

We left some open points, specially in section Indeed, it would be very nice to get a 
more complete understanding of Seiberg duality, the tension of the domain walls and the 
field theory living on them. Improving our understanding of the gauge theory, following 
for example the ideas in [21], would be desirable. It would be very nice to explore some 
hydrodynamical aspects of the solution dual to the field theory with Nf = 2Nc at finite tem- 
perature. Also, finding a black hole in the solutions with arbitrary Nf and N^ is complicated, 
but might prove very useful. It would be interesting to follow the formalism in [ZD] for our 
background, in order to construct an effective 5-d action with which the computation of any 
correlation function might be feasible (a lot of work should be done in this direction). Also, 
completing the approach we sketched in section [3 may have interesting consequences. 

Finally, it would be of great interest to find a sigma model describing strings in our 
backgrounds. In this line, the background with Nf = 2Nc, due to its simphcity, is the most 
suitable to start studying these issues. 

We believe that the solution to the problems presented above, among others, could signif- 
icantly improve the understanding of gauge-strings duality for QCD-like theories. 
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A Appendix: The (abelian) BPS Equations from a Su- 
perpotential 

It is instructive to obtain the BPS systems of section IHl using an alternative approach, 
computing the so-called superpotential. Obtaining the non-abelian system of section |3] with 
this technique would be much harder since it involves algebraic constraints. 

Generically, for a lagrangian of the form [i] is the "time" variable. A, (p°- a = 1,2, . . . are 
the fields): 



if one can find W such that: 



,dW dW 



„2 2 



2 dip'' d(p^ 4c2 

then the equations of motion are automatically satisfied by the first order system: 

dip'' , „„,dW 



drj 2c2 



(A.l) 



(A.2) 



(A.3) 



drj dip'' 

Now, let us recover the system ()3.19|) - ()3.22|) . The relevant action in Einstein frame is ()3.11|) . 
i.e. ()3.12j) + ()3.14|l . since ()3.15j) does not depend on the metric nor on the dilaton and therefore 
does not enter the Einstein equations (its effect was taken into account when writing ()3.18|) ). 

In order to obtain an effective lagrangian, we substitute the ansatz ()3.2|1 . ()3.18j) . We also 
impose by hand the condition = 4/. After integrating by parts the second derivatives 
arising from the Ricci scalar and defining 4Y = 8f + 2g + 2h + k, we find the effective 
lagrangian: 



sin 6 sin 



16Y 



12 



2h" - k 



12 



V 



(A.4) 
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with: 

8 8 ^ 

(A.5) 

The last term of this potential comes from the Dirac brane action (j3.14p after using (|3.16p . 
The lagrangian ()A.4|) is of the form ()A.1|) and, choosing C3 = 2, eq. ()A.2|) reads, 

V = lidgW)' + \{dhWf + {dkWy - (A.6) 
Remarkably, this equation has a very simple solution: 

W = -^{N^ - N})e-^^-^ + N^e-^^-^ - \^-^^+>' - g-^s+fc _ 2e~'' . (A.7) 

In particular, it is amusing to see how simple is the modification of this W when the flavor 
branes are added. Inserting ()A.7j) in ()A.3|1 (with the upper sign), we recover ()3.19|1 - ()3.22|) . 



B Appendix: The BPS Eqs. A Detailed Derivation 



In this appendix we analyze in detail the computations that led us to the system of eqs ()4.14j) - 
(|4.17p when flavoring the non-singular solution (|2.1|) - (j2.8p . The analysis of the generalized 
unflavored ansatz of section H?T] is incorporated in the same formalism just by taking Nf = 0. 

It is important to notice that the ansatz for the metric we will use in the following is a 
subcase of that considered in |21] so we can profit from their analysis (although Nf = in 
their case). For the sake of clarity, we rewrite here the ansatz of the main text ()4.9|1 . ()4.12|1 . 
We take the Einstein frame metric^^: 



dx'l 3 + dr' + e^^^'\de'^ + sin^Od^^) + 



2g(r) 2k(r) 

ed^f) + ^p(^3 + cosOd^f^ . (B.l) 

The vielbein we consider for this metric is the straightforward generalization of ()3.H1 by the 
inclusion of the a(r) dependence in and e^. Apart from the dilaton, there is also the RR 
3-form field strength: 

F(3) = -2Ar,e-3/-29-V Ae2Ae3 + ^6'e"=^/-3-VA(-e^Ae^ + e^Ae^) + 

+ ig-^Z-s^-'^ (Ar,(a2 _ 2ah + 1) - Nf) A A + 

+ e'^f-^^'-^N,{h - a) (-e^ A + A e^) A (B.2) 



^^Comparing to 24 (we write a subindex B for quantities defined in that paper): hig = /12B — Xb 
Kb — Qb = 0. Noticing that in the metric is written in string frame: Ab — 2/, dtB — 2e^^dr, 



ie4/+'i+5^ e9B ^ 2e''"f, e^P^ = %(,-S'S--^k-h-g ^ Afterwards, we wiU define g = and p = tij/2. 



2 



Notice also the definitions duB — e "^P^dtB, vb = e^P^+^^^ and that in "Jl, Nc has been set to 2. In order 
to use the most usual convention, we take a = —as, b = —Bb- 
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Let us analyze the dilatino and gravitino transformations. For the functions we have in the 
ansatz, they read: 
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The projections one has to impose on the Kilhng spinor are: 

e = ie* , Te^e = T^e , Trus^ = {A + B T^2)e 
Notice that consistency imphes: 



so one can parameterize: 
and then the spinor reads: 



A = cos a , ;B = sin a 



e = e 2 "^^eo 



where eo satisfies the unrotated projections: 



(B.3) 

(B.4) 

(B.5) 
(B.6) 
(B.7) 

(B.8) 



With this, it is straightforward to obtain the set of differential equations and constrains 
coming from ()B.3jl . The spin connection [de"" + a;"^ A = 0) is needed. We write it here for 
completeness: 



ip2 



fir 



4 
1 



= e~f (/' + g')e^ + -a'e ' ' "e^ 



23 



— e 



-S^k-2gA 



+ ae'f'^ cosh(A; - g)e^ 



UJ 



13 



~f+k-2g^2 



+ ae'f~^ cosh(A; - g)e'^ 



UJ 



tf>r 



.'^-f-h+g^v^ 



u =e -^{f + g)e - -ae '"e- , 
cj2i = (_e-/+'=-29 + 2e-f-^)e^ - e~^~^ cot ^e^ 



UJ 



3(p 



2h~f+k^^2 _ ^ ^^-f~h gij^j^^^ _ ^yl ^ 



,30 



UJ- = ^e-^''-f+\a^ - l)e^ + ae"^"'^ sinh(A; - g)e^ . 



(B.9) 



56 



Comparing the equation for the variation of the dilatino and that of Sip^, we immediately 
find the condition 

= 4/ (B.IO) 

Moreover, from St/jx — we get: 

/' = ^Ae-^<^-''-^Ae-^''-''{N^{a^-2ab+l)-Nf) + 

+^Be-''-''''N^{b - a) , (B.ll) 



2Nce^'-^-''B - l-e'-''-''B {N^ia^ - 2ab + 1) - Nf) 



-2Ae-^Nj^ - a) (B.12) 
From = (or equivalently, from Sip^^, = 0) and using the previous expressions: 

h' = -Be-'' (^a cosh(A; - g) + ^e-''-m,{b - a)^ - ^Ae-^''+\a^ - 1) + 

+^Ae-^''-'' {N,{a'' - 2ab + 1) - Nf) , (B.13) 

a' = -4e-Ma cosh(A; - g) + Be-^-'^+'=(a2 - 1) - 2N^Be''-^^-'' - 

—Be-''-^-'' (iVc(a2 - 2ab + 1) - Nf) (B.14) 

Prom 5-01 = (or 5V'2 = 0) : 

g' = ^e*^-'^ - Bae'^ sinh(A; - 5) - \Be-''-^-^Nj^b - a) - A^c^e-'^-'^ , 
= -4a^e'=-»-'^ - 4Be*^-2» + Be-2'^+'=(a2 - 1) (B.15) 

From (5-03 = one finds: 

A;' = ->le^-25 + 2^6"'= + i^e'2/.+^(a2 - 1) + 2ai3e-'' sinh(A; - g) - N^Ae'^^-k + 

+^^6-2'^-'= (Ar^(a' - 2ab + 1) - A^/) - Be'''-^-'^ N,{b - a) , (B.16) 
and the algebraic equation: 

= Be''-^<' - 2Be-'' - ^Be'^^^'' {a"^ - 1) + 2aAe-'' smh{k - g) + N^Be'^^''' - 

—Be-^''-'' {Nc{a^ - 2ab + 1) - Nf) - Ae'^'-^^-'^Ncib - a) (B.17) 
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Finally, 5ipr = yields: 

eo = e^r/, (B.18) 
where rj is some constant spinor satisfying the same projections as eo and: 

a' = - ^e-'^+^ (a' + N.b'e-^') (B. 19) 

In order to simplify these equations, let us make some definitions: 

and, following 

S = Ya^a^ + 2a2(-l + e29) + (1 + e^'^Y , 

C = — (l + a^ + e^s) . (B.21) 
2a 

Then, the set of algebraic constraints ()B.5jl . ()B.15jl . ()B.17|1 can be written as: 

A = — ^ , ^ = -^ (B-22) 



.29 



N,{bC - 1) + ^Nf 



aC-1 



(B.23) 



Although the system seems overdetermined, it is not. By deriving in ()B.22|) and using the 
rest of expressions, one arrives at ()B.19|) . which is therefore redundant. By deriving ()B.23|) 
and using all the differential equations and constraints, one finds an identity. Therefore, 
apart from ()B.19|) . one of the other differential equations (say the g' one) is also redundant. 
Let us count the number of functions. We have f,g,g,k,a,b. They are not independent 
because of ()B.23|) . But since the equation of g' plays no role, one can think of having five 
independent functions with a first order equation for each. 

Remarkably, the system of first order equations and algebraic constraints ()B.lH) - ()B.17j) 
solves the set of Einstein second order equations. The equations are ()3.23p and ()3.24|) along 
with the equation of motion for the RR form 9^ {^'^yJ—gioF^^^P^ = 0. We have checked this 
using Mathematica. 

We then would like to solve, after substituting the expressions (lR22|l . (lR23|) . the equations 
dEHH), (fRJ2|) . (fRl4|) . ^HJ^ and the equation for g that can be obtained from the above: 

g' = -Bae-^+3-^~Ae^^^+^{a^-l)-Ae^-^3+^Ae-^^-^ {Nc{a^ - 2ab + 1) - Nf)+N,Ae-^3-\ 

(B.24) 

Let us define: 

dp = e-^dr, (B.25) 
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then, one can prove that: 

dpC = 2S , dpS = 2C (B.26) 

and by also noticing the identity — S"^ = 1 and by fixing an integration constant (the 
origin of p) we can write: 

C = cosh(2p) , S = sinh(2p) (B.27) 
From this one can obtain for instance g in terms of a: 

e^^ = 2a cosh 2p - 1 - (B.28) 

The equation for b is reduced to: 



which is exphcitly solved by: 



2-^ 

^ N, 



7- j p + const 



sinh(2p) 



[B.30) 



We will fix this constant to zero. In the usual unfiavored case, setting the analogous constant 
to zero is required in order to have regularity at the origin. As we have seen in section |3J the 
fiavored case presents a singularity anyway. However, setting the constant to zero seems to 
allow a milder singular behavior at the origin. Since / never appears in the right hand sides 
of the differential equation, we are left with having to solve two coupled equations for a and 
k. These are the equations ()4.15|1 and ()4.16j) of the main text. It does not seem possible to 
find an explicit general solution of these equations. The equation for / is written in ()4.17|1 . 

Finally, let us comment on three particular explicit solutions of this system: When Nf = 0, 
one has, of course, the usual unfiavored solution of |Hl|7]: 

e- = l, a=^, e^/ = e^/o^ (iV, = 0) . (B.31) 

smh2p' 2^(pcoth2p--^-i)r ^ ^ ^ ^ ^ 

When Nc = Nf = one finds a Ricci fiat geometry which is nothing but the deformed 
conifold (see appendix for more details). And finally, when Nf = 2Nc, one has the 
solution presented in section 14.51 Notice that for these two last cases the manipulations 
leading to ()B.26|) - ()B.30|) and ()4.15|) - ()4.16|) are ill-defined, but the solutions can be obtained 
from the system (lRTTll - (lRT7j) . 



59 



C Appendix: Solutions Asymptoting to the Conifold 



When Nf = Nc = 0, the Ricci-flat solution compatible with the ansatz ()4.1|1 is the deformed 
conifold. Actually, (|4.H) is the natural way of writing the deformed conifold metric when 
it is obtained from gauged supergravity [7T]. The different functions of the ansatz take the 
following values: 

cosh(2p) 3 

= \el'^^^^^IC{p) , e'^ = el cosh(2p)/C(p) , (C.l) 
4 cosh(2p) 

where e is the deformation parameter, the dilaton is constant and we have defined: 

(sinh(4p) - 4py/' 
^^P^ = 2V3sinh(2p) • (^-^^ 

We have argued, and it is apparent from figure El that the unflavored equations ()4.6|) - ()4.8|) 
as well as the flavored equations ()4.14j) - ()4.17jl have large p solutions in which the effect of 
the branes becomes asymptotically negligible and they must therefore approach the Ricci 
flat geometry ()C.1|) . In this limit, the only difference between flavored and unflavored is a 
rescaling of e^^ , e^^, e^^ by (1 — 

In fact, this can be proved explicitly by discarding exponentially suppressed terms and 
considering an expansion starting with: 

acosh2p- 1 = /3e-^(2p- 1) + . . . (C.3) 

We can write: 

a = 2e-2^(l + /3e-^(2p-l)) + ... , = + . . . (C.4) 

and similar expressions for the rest of functions which agree at leading order with (jdlj) . By 
comparing ()C.4|) to ()C.1|) . we see that /3 is a constant related to the deformation parameter 
of the conifold. 



D Appendix: An Interesting Solution with Nf = 2Nc 

The differential equations ()4.15|) and ()4.16|) are not valid for the special value x = 2. To 
study the solutions in this case start then from the unsimplified equations ()B.1H1 - ()B.17|1 . 
In this appendix, we will show that there is another solution apart from the one presented 
in section 14.51 The manipulation of the differential equations and algebraic constraints 
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goes along similar lines to the x < 2 case, with the important exception that now we use 
the constraint ()B.17j) to solve for a and h rather than for e^^. We impose that the two 
independent terms in the equation vanish simultaneously, which gives 

^ and 6 = (D.l) 



cosh 2p 

The choice of this solution to ()B.17|) is suggested by the behavior of the x <2 solutions for 
X that goes towards 2, see the plots at the bottom of figure [T] 

At this point, it is easy to show that some of the remaining equations can be used to 
evaluate 

^ = tanh2p B = e''"^ = - tanh2p (D.2) 

^ cosh2p 2 ^ ^ ' 

and now, as before, we can reduce our system of BPS equations to a system of two coupled 
first order differential equations 

fD 31 

dpk = -{e^'' + Nc)e-''-3 + 2coth2p ^ 



with the only remaining unknown function determined once we find the solution to (jDT 

d,f = ^e-^-' (D.4) 

To solve equations ()D.3|) . we can obtain e^*^ from the first one, and substitute it into the 
second one to obtain a second order differential equation containing only e^~^^ 

^^9+hg2^g+h ^ ^Q^^g+hy + dpe3+''{3N^ - 2e3+^ coth 2p) + 2{N^- 6^+^ coth 2p) = (D.5) 

In the large p region, we can approximate coth2p ^ 1 up to exponentially suppressed terms, 
which gives us that e^"*"^ = Nc is an exact solution at large p. This also fixes the value of e^'^ 
in the large p region, e^'^ = Nc, and again the only corrections to this large p solution are 
exponentially suppressed. From ()D.2|) and ()D.4|) we find also 

~ 2Nc + 0{e-^P) (D.6) 
/ - 7 + • • • 



4 4v/iV, 

Notice that this behavior is very similar to the one of the other flavored x ^ 2 cases we 
described in subsection 14.3. even though the constants at infinity do not match. 

We can now plug these asymptotic conditions in Mathematica and obtain the numerical 
solution to the BPS equations. The functions we obtain are shown in figure El Notice that 
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p 



p 



0.5 1 1.5 2^ 0.5 1 1.5 2 2.5 3 



Figure 5: Behavior of the background functions for the x = 2 solution with a = ^^J^ . From 
the top left we have e'^^^p\ e^^^p\ e^3^p\ e^U{p)~h)^ Here we have taken = 1. 



this solution has a very nice property. In section 15.71 we had introduced the gravity dual 
description of Seiberg duality for the Af = 1 SQCD-like theories we are considering. The 
solution we have presented in this subsection is exactly Seiberg duality invariant. In fact, it 
is easy to check that this condition imposes a single equation 

e"" + ^{a" - 1) = (D.7) 
which is satisfied exactly by our a = 1/ cosh2p and e^~^ = |tanh2p. 

E Appendix: Some Technical Details in Section El 

In this appendix, we will just cover the technical details that were left aside in the analysis 
of gauge theory aspects of the unfiavored solutions. 

Massless Glueballs 

Let us start with the fluctuated eqs of motion, that lead us to the fact that a normalizable 
massless glueball does not exist. 
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The eqs for the fluctuation were analyzed in [58] and read, for the fluctuated Ricci tensor 
equation 



4 

^'"^ [g^,{2FsSF3 - S/i'^'F.^Fr + + h^.F^] (E.l) 
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For the fluctuation of the Ricci scalar, 



24 

For the fluctuated dilaton equation 



^'""^ [2F,5F, - 3h'^'F,,,Fr + ScjyF^l (E.2) 



^^S(l)-h^''V^d,(l)-^g"'dr(l)i2V''hr^-drh>;^)-^i6(l)F^ + 2F = (E.3) 

For the fluctuated Maxwell and Bianchi eqs. 



d^[^e^{{-hP + 5(t))F^'^^ + (JF'^^" - K'^F^" + K'^F^^ - /i^^F^")] = 0, d5F = 0. (E.4) 



,1 
'2 

One can easily see that a solution to these eqs is the one found by the authors of 

5Fs = *,da, h^, = 5(P = 0; Fs.5F; = (E.5) 

the analysis of this solution shows that is non-normalizable, hence is not a mode in the dual 
gauge theory. The rest is covered in the main part of the paper. 

Rotating strings 

Now, let us turn to the analysis of new solutions for strings rotating in our backgrounds. 
So, the Polyakov action for the configuration we will propose reads (r, a are the worldsheet 
coordinates with the usual metric diag(-l,l)) reads, 

5 = ^ y dadrG^^daX^dbX^r]'''' (E.6) 

and the Virasoro constraints: 

drX^d^X^ = 

G^.drX^drX" + G^.d^Xf^d^X" = (E.7) 
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Consider the ansatz 



t = KT, X = x{T,a), y = y{T,a), 

p = p{a) , ^ = V>(r) (E.8) 
This is consistent ii 6 = 6 = ^. We consider the induced string frame metric 



ds"" = i^-de + dx^ + dy^ + e^'^dp'' + ^d^j^ j (E.9) 

where and k{p). The equations of motion are consistent with: 

X = -R(cr) cos(ti;ir) , y = R{a) sin(a;ir) , ip = 2uj2T (E.IO) 

One obtains (prime is derivative with respect to a): 

R" + {dp(j))p'R' + ujIR = 
p" + {d,k)p'^ - (9,0)e-2'=i?'2 + {d,k)ul = 

-K^ + ulR^ + R!"^ + e^^ul + e^V^ = (E.ll) 

These equations are not independent since by deriving the last one and using the other two, 
one reaches an identity. The rest follows as in section |H1 

Penrose Limit and PP-waves 

Now, let us present the technical details of the Penrose limit and plane wave geometry. 
The computation is very similar to the one in [U^ . 

After doing a gauge transformation (that can be better understood as a coordinate trans- 
formation or reparametrization of the left invariant forms of SU{2), see we leave the 
geometry in the form (j4.H) . but with the S'f/(2)-valued one form y4"cr" being 

A = pp^ (cos LpdO — cos 6 sin 6 sin Lpd(f)a'^ + (sin LpdO + cos 9 sin 9 cos <^dLp)a'^ + sin^ 9dLpa^ . 

(E.12) 

Now, we have to rescale and redefine the coordinates according to, 

P = j. = j, f=|, 20+ = V + <^ (E.13) 

Then we take the limit of L — »• oo and keep only terms suppressed in L~^. We get a PP-wave 
that contains the so called magnetic terms (that is crossed terms d(j)+d(p and dip+dip ), that 
reads, 

g-2/c(O)-0(O)/2^2^^2 ^ ^2^_^^2g-2fc(0) ^ ^^2^) ^ ^^2 ^ ^2^^^2 ^ ^^^2 ^^^2) ^ 

1 v"^ 

-{dv^ + v^dif'^) - {—d^ + pr"^ sm^ 9dip)d<i)+. (E.14) 
4 2 



,2fc 
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We remind the reader that /i is the parameter labelhng the different solutions near p = 0. 
Even when this can be quantized, is better in order to have a better intuition on the spectrum, 
to make a last change in variables 

dip dip + ^d(f)^, dip —>■ dip + d(j)^, (E.15) 
then, defining as usual, = te~^^^\ x~ = L'^{t — 0+), we have that the geometry looks, 

g-2fc(0)-<^(0)/2^2^^2 ^ _2dx+dx- + dxl + dyi + dz' + \dvl - + ^yl)dcfi- (E.16) 
For the rest of the analysis, please go back to the relevant part in section |H| 
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